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SECTION I 


INTRODUCTION AND SUMMARY 
1.1 INTRODUCTION 

The development of laser communication systems for space is 
a vital step toward providing the nation with the capability for synchron- 
ous relays, deep space probes, and ground stations. The special 
properties of laser systems make such links possible without the need 
for large antenna structures, and the need for high power as required 
by microwave systems. Such a capability is particularly required, 
within this decade, to furnish wideband data links from low altitude 
satellites with high resolution sensors to earth synchronous communica- 
tion terminals from which the information may then be relayed to 
ground stations. 

Spaceborne communication systems must necessarily operate 
near theoretical thresholds. The demands against weight, prime 
power, information rate, cost, reliability, and other factors do not 
allow a large safety margin to compensate for overlooked degradations. 
For this reason and others, the present computer program, LACOMA 
(LA ser C OM munic ato r Analysis program), has been developed to 
accurately simulate optical configurations of heterodyne communication 
systems and to evaluate the effects of aberrations, distortions and mis- 
alignments on the heterodyne performance of these systems. The pro- 
gram has been developed with the communication engineer in mind. It 
has been designed to be easy to program and easy to interpret. Output 
is specified in terms of antenna gain to allow results to be directly 
applied to system calculations. Absolute power values are also printed 
out to provide actual signal levels for evaluation. 
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1.2 


COMPARISON WITH OTHER OPTICAL COMPUTER PROGRAMS 


LACOMA represents a major advance in the optical analysis 
of laser heterodyne communication systems. It also represents an 
advance in general optical system analysis since it includes the effects 
of gaussian pupil functions. Some of the optical analysis programs with 
which LACOMA might be compared are: 

1. POLY PACOS — Aerospace Corporation 

2. SLAP (Spectral Lens Analysis Program) — PAGOS 
Corporation 

3. MTF— D. Grey Associates 

4. ACCOS-GOALS— Scientific Computations Inc. 

These represent a cross section of the generally available opti- 
cal analysis progrcims and will serve as benchmarks to compare 
LACOMA. 

LACOMA is a direct descendant of SLAP which is in turn an 
updated version of POLY PAGOS, so that one would logically expect 
more similarities between these programs. An overview of the pro- 
grams is provided by the tabular check sheet of Table It 1 . A major 
difference between these programs is the Fourier transform algorithms. 
LACOMA (also SLAP and POLYPAGOS) utilizes an algorithm which per- 
mits the specification of number and location of the output points for the 
computed spread function (ASF, PSF) results. For a given optical 
system or telescope, LACOMA allows sampling of the Airy disc via 
grid sizes from 2x2 to 101x101 or any other arrays between these 
extremes. This is not possible with the algorithms used by the other 
programs. Not included in the table is a cost comparison. If it were 
possible to set up a cost or running time comparison, the Gray-MTF 
might show somewhat better due to its heavy emphasis on machine 
language programming for the CDC 6600. The ACCOS-GOALS would 
probably score as the most costly. Some of the tabulated categories 
are subjective such as the ease of input, but they reflect more than a 
single view. 
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TABLE 1-1 


Comparison of Optical Computer Programs 


Feature 

LACOMA 

SLAP 

POLY- 

PAGOS 

MTF- 

GREY 

ACCOS/ 

GOALS 

Heterodyne Power 

X 





Detector Power 

X 

X 

X 

X 


Multiple Detectors 

X 





Multiple Beams 

X 





Optimization Capability 





X 

Interactive Operation 





X 

Gaussian Input Pupil 

X 





Transmittance 


r 




Output Grid Selection 

X 

X 

X 



Deformed Surfaces 

X 

X 

X 



Perfect Surface 

X 

X 

X 



3-D Tilts and Decentrations 

X 

X 

X 

X 

X 

Geometrical Analysis 


X 

X 


X 

Generalized Aperture Description 

X 

X 




Ease of Input 

X 

X 

X 


X 


T1199 
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LiACOMA constitutes a new tool in the optical communications 
field. It is unique in its ability to perform an accurate analysis of the 
laser communication system. 

1.3 RECOMMENDED ADDITIONS FOR INCREASED VERSATILITY 

LACOMA has been shown to be a unique and useful tool in the 
analysis of heterodyne systems (see Section III). During the course of 
using the program, it was often tempting to use it as a design tool. 

This led to an investigation into the possibility of modifying the present 
program to combine the features of design and analysis. The result 
of this study led to the following list of recommended additions, some 
of which would provide some design or optimization capability, and 
others that would enhance the analysis function. 

a. Optimization — Conventional optimization techniques will 
not necessarily produce optimum optics for laser communication sys- 
tems, especially for the transmitter and local oscillator. A gaussian 
beam amounts to an aperture weighting function which should be included 
in the optimization process. Also, ideally, the received signal optics 
and local oscillator optics should be optimized on the basis of their 
combined beams. It is recommended that work in the optimization 
area occur in stages. 

(1) Focus optimization — optimization of best focus posi- 
tion for a given system configuration. 

(2) General optimization Phase I — Development of 
algorithms and approach to optimization for trans- 
mitter and receiver. 

(3) General optimization — Phase II— Development of 
optimization program based on algorithms and 
approach spelled out in Phase I. 

b. Tolerance Analysis — Automatic determination of toler- 
ance budgets necessary to maintain specified level of performance. 

This can be a separately addressable subroutine of an optimization 
program. 
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c. Perturbed Focal Surface — Revisions to permit tilted, 
decentered, curved focal surface without dummy surfaces. 

d. Return Option for Tilted- Dec entered Surfaces — To auto- 
matically restore coordinate system after tilt and/or decentration 
without dxammy surfaces. 

e. Aberrated Entrance Pupil — Accommodation of aberra- 
tions of entrance pupil due to tilts- dec entrations, input beam inclination. 

f. Potato Chip Aspheres — Generalize deformed surface 
options to include more generally deformed surfaces due to thermal- 
mechanical stresses or potato chip aspheres. 

g. Pseudo- diffraction Intensity Function — Geometrical 
analysis to be used for evaluation of systems without involving full scale 
diffraction- heterodyne analysis. 

h. Transmittance — Subroutine to compute losses due to 
absorption losses, Fresnel reflection losses, antireflection films, 
multilayer thin films. 

Polarization — extension to determine orthogonal polarization 
components of beam. 

i. Graphics — Isometric and contour representation of wave 
fvinction, pupil function, spread functions, etc. 

Interactive graphics — for analysis and optimization. 

j. Fresnel Diffraction — Modification of diffraction integrals 
of LACOMA to permit computation of amplitude or intensity function out 
of Fraunhofer plane. 

k. Thermal-mechanical Structural Analysis — Extension of 
LACOMA to interface with thermal-mechanical analysis programs such 
as STARDYNE, NASTRAN, etc., to perform overall system analysis. 

Some of the listed modifications are contained in other programs 
as seen in Table 2-1, while others, such as transmittance, are unique. 
Transmittance in particular shotild be considered for LACOMA, since 
this parameter affects the amplitude distribution across the exit pupil, 
which, in turn, can make a significant difference in phase distribution 
across the focal plane, and can also cause a shift in the location of the 
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focal plane. Other modifications that are considered particularly 
appropriate to heterodyne systems are focus optimization, perturbed 
focal surfaces, and the related Fresnel diffraction modification. 

1.4 SUMMARY 

Section II of this report is designed to be a self-contained user's 
manual. An engineer familiar with heterodyne optical systems should 
be able to utilize LACOMA with this manual and a minimal amount of 
help from computer personnel. 

Section III contains an analysis of the sample system, the opti- 
cal mechanical subsystem developed for NASA under contract NAS- 
5-21859. This analysis includes a "perfect" reference system, several 
runs for the system with no distortions or errors, and finally analyses 
with tilts and distortions that might be introduced by manxifacturing 
tolerances or by mechanical or thermal stress. This section has been 
arranged to provide a specific example for using LACOMA. 

At the end of Section III, two additional analyses are summarized, 
one that establishes bounds on the field of view of heterodyne receivers, 
and a second that analyzes a telescope designed for a terrestrial system. 
This last example was selected because there is sufficient aberration to 
cause significant performance degradation, thus providing a means for 
illustrating the advantages of LACOMA over conventional optical 
programs. 

Section IV contains a detailed discussion of the program. 

Specific equations and computational flow diagrams are discussed. 

Error sources and their probable magnitudes are covered at the end 
of this section. 
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SECTION II 


LACOMA USER'S MANUAL 

2. 1 GENERAL PROGRAM DESCRIPTION 

The program has been developed to allow prediction of laser 
communication system performance and to permit tradeoff analysis in 
the design of these systems. The program is set up to allow the 
optional analysis of transmitter or receiver optics. 

For the transmitter, the program starts with the data for the 
specified laser beam and propagates the beam through the optical train 
to determine the far-field intensity function. 

For the receiver, the program carries the received signal beam 
through the optical train to the detector. The local oscillator laser 
beam is also traced to the detector where the two beams are combined 
and the various quality criteria computed. 

The computation for the receiver requires the entry of two sets 
of optical data — one for the received signal optics and another for the 
local oscillator optics. The program checks whether one or two sets 
of data have been entered (see Fig. 2-1). If there is only one set, it 
performs the transmitter analysis. Two sets of data precipitate the 
receiver analysis. 

a. Transmitter — For a single set of optical data the pro- 

gram computes the gaussian beam profile and any specified quadratic 
phase on the input beam and traces it through the optical train, carrying 
the complex cimplitude function through in the form of an amplitude 
array B (i, j), and the optical path difference array W (i, j) (see Fig. 2-2). 
At the exit pupil it computes the rms wavefront error crW and the 
normalized Strehl intensity I (cr) for the (systematic) system wavefront. 
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DATA INPUT 



Fig. Z-1. Basic flow diagram. 
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Fig. 2-1 
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(^DENOTES FOURIER TRANSFORM) 


Fig. 2-2. Transmitter computation flow 


NAS 5-21898 

Final Technical Report 

Fig. 2-2 
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It forms the pupil function F (u, v) (= complex amplitude array at the 
exit pupil). 


F(u, v) = B(u, v) • exp [-KW(u, v)] 



u, V = pupil coordinates 


(2-1) 


The Fourier transform of F(u, v) yields the far-field amplitude 
spread function ASF(x, y) 


ASF(x,y) = ^ F(u, v) (2-2) 

^ denotes Fourier transform. 

Squaring the modulus of the complex amplitude spread function 
gives the far-field intensity or point spread function PSF (x, y) 

PSF(x,y) = |ASF(x,y)|^ (2-3) 

b. Transmitter Quality Criteria — The far-field intensity 

function is itself a measure of the transmitter performance. This is 
augmented by the computation and output of: the peak far-field intensity 
value, optical transmission, maximum antenna gain, and the overall 
transmitter efficiency. 

• Peak intensity — This is simply the peak value of the 
far-field intensity function and is included in the 
quality criteria in case the investigator chooses 

not to print out the entire far-field intensity 
function. 

• Optical Transmission — This is the ratio of the 
exit pupil power P^ to the total laser power P^ 

P 

Optical Transmission = (2-4) 

o 
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Transmission losses accovinted for here are those due 
to vignetting and obscuration. 

• Maximtim Antenna Gain — This is computed for the 

exit pupil area A^. 


4-rTA 

0 

Maximum Antenna Gain = r— 


(2-5) 


• Overall Transmitter Efficiency — This value relates 

the computed peak far-field intensity of the trans- 
mitter to that for an ideal transmitter with the same 
input power and effective exit pupil area 


Overall Transmitter Efficiency = 


2 2 
I R X 
r 

A P 
e e 


(2-6) 


= Transmitter peak intensity 
R = Range 

c. Receiver — When two sets of optical data are entered, 

the program performs the receiver analysis. The optical data are 
entered sequentially. That is, there is a total of N surfaces entered 
where 


N = Nj + N2 

The first set of surfaces represents the received signal 
optical train and the following surfaces are for the local oscillator 
optical train. The progrcim (Fig. 2-3) computes the received signal 
pupil fvinction F^^(u, v). 


11 


3079-3 



Fig. 2-3. Receiver computational flow. 
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Fig. 2-3 
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( 2 - 12 ) 


Nonheterodyne Detection Efficiency 



P 

o 


• L.O. illumination Efficiency_— This efficiency com- 

pares the heterodyne power P, for the two given 
signals with no phase error with the optimum 
heterodyne power value for two signals of cor- 
responding power 


Li. O. Illumination Efficiency 




(2-13) 


• Maximum Antenna Gain — This is the theoretical 

gain for the effective entrance pupil area A of 
the received signal optical train. 


Maximum Antenna Gain 


4ttA 


(2-14) 


• Receiver Efficiency to i.f. — Assuming a quantum 

. efficiency of 1, this efficiency relates the square 
of the heterodyne power P^ to the product of the 
total received power Pq and the L.O. power Pl, 
and is a measure of overall receiver performance. 


Receiver Efficiency to i.f. 


P P, 
o L 


(2-15) 


2 . 2 DATA DISC USSION 

The program performs analyses on laser communication system 
optics as summarized in Section 2-1 with the detailed computational 
procedure described in Section IV of this document. The analyses per- 
formed by the program are briefly listed here. 

Paraxial Analysis (Section 4. 1) 

Ray Trace-Optical Path Difference (Section 4.3) 
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Amplitude and Point Spread Function (Section 4. 8) 

Receiver Quality Criteria. (Section 4.9) 

Transmitter Quality Criteria (Section 4. 10) 

Numerous options are available for handling of the data, choice 
of computations, output data, etc. Two data arrays IPR(?G and IDETEC 
control the major f\inctions and are discussed in detail. 

Some optical system parameters and/or configurations which 
can be accommodated by LACOMA include: 

• Catoptric systems 

• Catadioptric systems 

• Dioptric systems 

• Spherical surfaces 

• Aspheres 

• Cylindrical surfaces 

• Toroidal surfaces 

• Surfaces with slight cylindrical error or warpage 

• Periodic surface errors 

• Misaligned (tilted, and/or decentered) systems 

• Composite vignetting and obscuration 

• Afocal systems (via perfect imaging lens) 

Any of the analyses can be performed for any of the system 
parameters or configurations. 

a. Data Requirements — A certain minimum amount of data 
must be entered to obtain any results. These minimum data require- 
ments can combine with built-in data or default conditions to give a 
complete system analysis. These data are listed here with their 
associated program variable names and text references. 
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Data 

Niimber of surfaces 
Surface curvature 

Surface- surface spacing 

Refractive index 

Entrance pupil semi- diameter 


Variable Name 
N 

RHOS (p ) 
s 


TS (t ) 

' s 


XMUS (Pg) 
BET AO ((3^) 


Text Reference 

2.6, Fig. 2-1 

2.6, 4. 1, 4.3, 

Fig. 2-1, 4-1 

2.6, 4. 1, 4.2, 4.3, 
Fig. 2-1, 4-2 

2. 6, 4. 1, 4.3, 

Fig. 2-1 

2.6, 4. 1, Fig. 2-1 


With these data input, quality criteria will be obtained for a transmitter 
or receiver depending on the entry of values for N(l) only or N(l) and 
N(2). The computations will be performed under the following default 
conditions. 


(1) Plane wavefront input- Uniphase energy 

(2) Analysis performed at paraxial image position 

(3) Axial analysis only — no inclination of input beam 

(4) Entrance pupil at first surface 

(5) No misalignment 

(6) No vignetting or obscuration 

(7) All surfaces rotationally symmetric spheres or aspheres 

(8) Wavelength = 0.010611385 (mm) 

(9) Data ass\xmed input in millimeters 

(10) For receiver — single circular detector, centered 
on optical axis, detector radius = radius of AIRY 
disc for uniformly illiiminated unobscured aperture 
of specified semi-diameter. 

(11) Pupil array = 51 x 51 

(12) Spread function array =51x51 

(13) Detector array = 21 x 21 (round) 

(14) Output data = Paraxial analysis, wavefront statistics, 
quality criteria 

2 

(15) Gaussian beam trvincation point assumed l/e intensity 
point. 
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b. Data Discussion — System parameters and program 

options effected by the default conditions are listed below. The rele- 
vant variable names available to override the defaxolt conditions are 
given in parentheses with page references for descriptions of the 
variables. 

(1) Input wavefront or object distance Section 4. 1 , 4.7 
XLINV (Lo'^) 

MAG 

(2) Image (or detector) location Section 4. 2 

IFLGl 

TS(N-l) (tj^_j) 

TS (N) (t^) 

(3) Inclination of input beam Section 4.7 

H 

(4) Entrance Pupil location Section 4. 1 

M 

(5) Misalignment Section 4. 6 

DLTAX (^x) 

DLTAY (Ay) 

DLTAZ (Az) 

THTAX (0x) 

THTAY (0y) 

THTAZ (0z) 

(6) Vignetting or obscuration Section 4. 4 

BETASX (6 ) 

''^sx 

BETASY (B ) 

'•^sy 
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BETAPX (B ) 

BETAPY (B ) 

^py 

VZER0 (V^) 

VPI(V^) 

GZER0 (G^) 

GPI (G^) 

UZER0 (U ) 
o 

(7) Special surfaces Section 4. 3 

RH01 (Pj) 

RH<Z2 (p^) 

CC (b) 

CRS (a) 

PHI (4>) 

AMP (C j ) 

FREQ (C^) 

C3 (C 3 ) 

USI (U^) 

VSI (V 3 ) 

( 8 ) Wavelength Section 4. 8 

VLAMDA (\) 

(9) System units of length^ All Dimensional Data Must Be 
Input In Same Units as Wavelength. 

(10) Detector Parameters Section 4 . 8 

IDETEC 

NDET 
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( 11 ) 


Section 4. 7 


Pupil Array 
NXY 

(12) Spread function array Section 4.7 

. NFRS 

DELFS 

(13) Detector Array 

IDETEC Section 4. 8 

DELFS 

(14) Output Data 
IPR0G 

(15) Gaussian beam parameters Section 4.7 

Z0MEGA (oj^) 

PZER0(P^) 

2.3 SETTING UP THE DATA DECK 


The data deck consists of the following: 

• Title Card — one card with any combination of 
characters in columns 1 to 80, 

• First Data Card — The first card must contain ^INP 
in columns 2 to 5, followed by a blank. Data may 
begin on this card or the next. 

• Data Cards — As many as necessary with data in 
cols. 2 to 80. 

• End of data card — the last card must contain 
$END in cols. 2 to 5. 

The details of the input format of these cards are discussed in 
the Input Format, Section 2.5. 
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The details of the individual data are given in Input Parameters, 
Section 2.6. 

2.4 DATA PREPARATION PROCEDURE 

For the receiver analysis, two sets of data are entered sequen- 
tially. The first N(l) surfaces are the received signal optics. These 
are immediately followed by N(2) surfaces representing the local 
oscillator optics. A total of N(l) + N(2) surfaces are entered, surface 
No. 1 is the first surface of the received signal optics, surface No. N(l) 
is the last (detector) surface for the received signal optics, surface 
No. (N(l) + 1)) is the first surface of the local oscillator optics and 
surface No. (N(l) + N(2)) is the last (detector) surface for the local 
oscillator optics. The system parameters for the received signal 
optics versus the LO optics are defined by two component arrays with 
the first component representing the received signal optics in all cases. 

Examples: The sample case entitled "RECEIVER TEST CASE 
#2", (Appendix D) indicates surface data for a total of twenty 
surfaces. Also input are 

• N(l) = -12, -8, where the negative signs indicate 
that data is input in radius form rather than 
curvature, the first twelve surfaces are the 
received signal optics and the following eight are 
the LO optics, thus accounting for the twenty 
surfaces. 

• BETAO(l) = 53. 975, 2.621, the received signal 
optics semi-diameter is 53. 975 mm and that for 
the LO is 2.621 mm. 

Notice that, as discussed in Input Format, Section 2. 5, 

N(l) = -12, -8, 
is the same as 

N(l) = -12, N(2) = -8, 
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and 


BETAO(l) = 53.975, 2.621, 
is the same as 

BETAO(l) = 53.975, BETAO(2) = 2.621, 

The program indication of whether a transmitter or receiver 
is being analyzed is based on the value of N(2) 

If N(2) = 0 — Transmitter analysis 
If N(2) 9 ^ 0 — Receiver analysis. 

It is recommended that the user refer to the facsimile input 
sheet, Fig. 2-4. All of the variables are shown on this sheet. 
The user can look at each variable, decide whether it applies 
to his problem, enter an appropriate value if it does and move 
on to the next variable. This procedure guarantees that no 
data is overlooked. Notice that the variables may occur in 
• any order in the data deck which may or may not correspond 
to the order in which the variables appear on the facsimile 
input sheet. 

2.5 INPUT FORMAT 

Input data consists of two groups, both of which must be present 
for each case, in the following order: 

a. One title card — any combination of blanks, letters, 
numerals, and the characters + -, =/$(). * in 
colvimns 1 to 60. 

b. One or more cards of nvimerical data in colums 2 to 80, 
which are read in by means of the NAMELIST input 
feature, as described below. 
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$INP N(1) = , , M(l) = , , BETAO(l) = , , XLINV(l) = , , MAG(l) = , , 

IFLGI(l) = , , H(l) = , , RHOSO) = , , . , TS(1) = , . , . XMUS(l) = , , , 

SKAPA(l) = , , , DFORM(l,l) = , . , RADS(1) = , , , RH01(1) = . , , 

RH02(1) = , , , PHI(l) = , , , AMP(l) = , , , FREQ(I) = , , , USI(l) = , , , 
VSIO) = , , , C3(l) = , , , CRS(l) = , , CC(1) = , , THTAX(l) = , , , 
THTAY(l) = , , THTAZO) = , , DLTAX(l) = . , DLTAY(l) = . , OLTAZ(l) = , , 
BETAPX(I) = , , ,.BETAPY{1) , , , BETASX(l) = . . , BETASY(l) = , , , 

VZERO(l) = , , , VPI(l) = , , GZERO(l) = , , , GPI(1) = . , UZERO(l) = , 

VLAMDA = , NFRS = , DELFS = , NXY = , PZERO(l) = , , ZOMEGA = , 

IPROG(l) = , , , IDETEC(l.l) = , , , NDET = . FLAG = , 

$END 


Fig. 2-4. Facsimile input sheet. 
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NAMELIST REQUIREMENTS 

a. Column 1 must be blank on all cards, 

b. The first card of the data must contain $INP 
in columns 2 to 5, followed by a blank space. 

Variable values may begin on this card or on 
the next card. The termination of a data set 

is indicated by the characters $END in columns 
2 to 5. The appearance of the character $ any- 
where else will cause an error. 

c. Variable names 

1. Single -valued variables are input in the form 

(name.) = (value) 

Example: NPTS = 16, DX = 0.02, 

2. Commas follow every numerical value. 

3. Arrays may be put in as single elements with 
the subscript or by listing the consecutive 
values: 

Examples: A(l) = 3, A(4) =-1.69, A(5) = 32.1, 
or 

A = 3, 0., 0., -1.69, 32,1, 

4. Double- subscripted arrays must be input 
columnwise. 

Example: A 4 x 3 matrix A may be input in the 
order 

~ ^ir ^21’ ^31’ ^41’ ^12’ ^22’ ^32’ ^42’ ^13’ 
^23’ ^33’ ^43 

5. A string of consecutive elements in an array may 
be entered by giving the name and subscript of 
the first element. 

Example: To input elements 19 through 23 of 
array A, write 

A(19) = ajg, a^Q, a.2i> ^22‘ ^23 
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d. Value formats 

1. Whole numbers may be input with or without a 
decimal point. Exponents (power of ten) may be 
indicated by an E followed by the power, or the 
E may be omitted and a signed integer used 

for the power. 

7.46 X 10^ may be written 7.46E6, 
7.46E+6, 7.46+6, or 
7.46+06: 

7.46 X 10'^ may be written 7.46E-6, 
7.46-6, or 7.46-06 

2. No plus signs are necessary for positive values 
or exponents. Negative values or exponents are 
indicated with a minus sign. 

- 8 

Example: -4.396 x 10 becomes -4.396E-8 
or -4. 396-8 

3. Double precision numbers take a D instead of 
an E to indicate the exponent. If a double 
precision number contains less than nine sig- 
nificant digits, it must have a D plus exponent 
for proper conversion. 

Examples: 3. 141592653587973 is written just 
like that, while 3. 14159 becomes 
3.14159D0. 

4. Identical consecutive values of an array may 
be abbreviated by writing an integral multiple 
and an asterisk(^') in front of the value. 

Excimple: if = 2, and A^ = A^ . . . 



write 

A = 1., 2., 31*-4, 
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e. Errors 

1. If the $INP followed by a blank does not appear 
as the first punches in the first card (excluding 
column 1), the computer will ignore that card 
and continue reading, trying to find $INP in the 
next cards. This process continues until $INP 
is found or there are no more cards. No error 
message is given if the wrong cards are being 
read and rejected. 

2. If a variable name is misspelled, the computer 
will give an error trace, terminate execution, 
print the following message; "Namelist name 
not found. " 

f. Order multiple cases 

1. Variables may appear in any order. 

2. Not all data need appear in any set. On suc- 
cessive cases where it is desired to change 
just a few values, only those variables need 
be input, with the rest retaining their values 
from previous cases. 

3. Values are all zero initially, unless other- 
wise specified. Thus, it is necessary to 
input only nonzero values. If runs are 
stacked, however, any data not written over 
will carry over from the preceding run. 

4. The same variable name may appear two or 
more times in a data set. The value physically 
last in the deck will override any previous 
values. Thus, it is not necessary to repunch 
cards to change n\imbers, just place a card 
with the changed value (and variable name 
with proper subscript, if any) somewhere 
following the old value. 
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2.6 INPUT PARAMETERS 


Variable 

Name 


N(I) 


M(I) 


BETAO(I) 


XLINV(I) 


MAG (I) 


Description, Comments, Text References 

One- dimensional array of length (2), specifying number 
of surfaces. Count the image surface but not the 
entrance pupil or object. A negative sign affixed to 
N indicates input of RADS (radius) instead of RH0S 
(curvature). If N(2) = 0, transmitter analysis is to be 
performed. If N(2) 0, N(l) specifies the number of 

surfaces in the received signal optics and N(2) is the 
number of surfaces in the local oscillator optics. 

One- dimensional array of length (2) for surface number 
of aperture stop. M(l) = surface number of received 
signal (or transmitter) aperture stop. M(2) = LO 
aperture stop surface number. Default values; 

M(l) = 1, N(2) + 1. 

One-dimensional array of length (2) for semi-diameter 
of entrance pupil. BETAO(l) — received signal or 
transmitter. BETA0(2) — LO. 

One-dimensional array of length (2) specifying inverse 
of radius of incident wavefront at entrance pupil. 

XLINV(l) — Received signal or transmitter. 

XLINV(2) — L.O. 

Default values: XLINV(l) = 0. , 0. , 

One- dimensional array of length (2) specifying mag- 
nification of output versus input wavefront radius. 

MAG - -Rn/R„ 

where 

Rj^ = output wavefront radius 
R = input wavefront radius 
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and the sign conventions on R^, R^^ are standard. 
MAG(l) — Received signal or transmitter. 

MAG(2) ^ L.O. 

Default: MAG(l) = 0., 0., 

IFLG1(I) One- dimensional array of length (2) controlling location 

of image (or detector). 

IFIFLGl = 0 T' - = BF + 

N- 1 N 

IFIFLGl = 1 + Tj, 

where 

= spacing used in analysis as spacing 
from surface (N-1) to image 

BF = Paraxial back focus 

j = input value TS(N-l) 

= input value TS(N) 

IFLGl(l) — Received signal or transmitter 
IFLG1(2) ^ L.O. 

Default: IFLGl(l) = 0,0, 

NOTE: For IFLGl = 1 and 0, j will not 

revert to its previous value when running 
consecutive cases. 

H(I) One- dimensional array of length (2) specifying obliquity 

in radians of input beam. 

H(l) Received signal or transmitter 

H(2) L.O. 

Default: H(l) = 0., 0., 
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RH0S(I) 


RADS (I) 


TS(I) 


XMUS(I) 


DF0RM(J,I) 


One -dimensional array of length (101) for the spherical 

f h 

curvature (or base curvature of an asphere) for the D 
surface when N is positive. The sign convention to be 
applied to RH0S is that the curvature of a surface is 
positive when the center of curvature lies to the right 
of the surface. 

One-dimensional array of length (101) for the spherical 
radius (or base radius of an asphere) for the optical 
surface. The sign convention for RADS is the same as 
that for RH<2S. 

N must be entered negative when data is input into this 
array. 

RADS(I) = 0 defines a plane surface. 

One- dimensional array (101) for the axial separation 
between surface I and I + 1. 

The sign convention for TS is that TS is positive for 
spacings measured from left to right. 

One- dimensional array (101) for the index or refrac- 
tion of the medium between surfaces I and I + 1. 

A sign change between XMUS (I- 1) and XMUS(I) indicates 
a reflection at surface I. 

Default: XMUS(I) = 1. 

Two-dimensional array (5, 101) for aspheric deformation 

coefficients of the 2(J+1) power terms at surface I 

DFORM(l,I) = a 

DFORM(2,I) H p 

DFORM(3,I) = y 

DFORM(4,I) S J 

DFORM(5,I) = ~ 
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The aspheric expression is 


SKAPA(I) 


RH0ia) 


RH02(I) 


c = 


PQ 


1 + (1 -K p^Q^)2 


+ O'Q'^ + (3Q^ + YQ® + 


where 



x,y = surface intercept points 
p = RH0S = base curvature 
K = (1 - € ) = conic coefficient — SKAPA 

One- dimensional array (101) for conic coefficient of 
surface I 


where 

SKAPA = K = 1 - € ^ 


£ = conic eccentricity 


£ 



+ S 


1 

1 


Sj = distance from first focus to conic 

= distance from conic to second focus 

One- dimensional array (101) for base curvature of 
aspheric (acircular) cylinder. Input in addition to 
RH0S(I). RH01 is used for tracing exact rays through 
surface I while RH0S is used in the paraxial ray trace, 
determining the first order parameters of surface I. 


One- dimensional array (101). Toric rotation curvature 
of surface I. RH02 used (with RH01) for exact ray 
trace. RH0S used in paraxial analysis. The optical 
axis and the axes of RH01 and RH02 are mutually 
orthogonal. 
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CRS(I) 

CC(I) 

pm (I) 

AMP (I) 
FREQ (I) 


C3a) 


One- dimensional array (101). Constant of rotational 
symmetry at surface I. CRS(I) = 1., for rotationally 
symmetrical surface of curvature RHQ1(I). CRS(I) 

= -1. identifies surface I as a "perfect" surface. 

Cylindrical constant. When CC(I) = 0 and CRS(I) = 0. 
surface is cylinder with curvature RH01 or RHQZ 
depending on which is nonzero. If RHQI and RHQZ 
are nonzero, surface is toric. 

One- dimensional array (101). Angular rotation of 
cylindrical axes. For Pm(I)=0 the axis of the RH01 
cylinder is parallel to the y axis and the axis of the 
RHOZ cylinder is parallel to the x axis. 

One-dimensional array (101). Amplitude term for 
periodic surface error function at surface I. 

One- dimensional array (101). Frequency component 
associated with sine function portion of periodic sur- 
face error. A surface of semidiameter BETAO 
will exhibit NZ siiic function zeroes between the center 
and the edge for 


NZ ^ 
BETAO ■ 


FREQ < 


NZ + 1 
BETAO 


One- dimensional array (101). Frequency component 
associated with cosine function portion of periodic 
surface error. A surface of semidiameter BETAO 
will exhibit NO complete cycles between the center 
and the edge for 

INC-. , 2JNC+12 

BETAO - BETAO 
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USI(I) 

VSI(I) 

DLTAX(I) 

DLTAY(I) 

DLTAZ(I) 

THTAX(I) 

THTAY(I) 

THTAZ(I) 

BETAPX(I) 

BETAPY(I) 


One- dimensional array (101). Displacement in 
X direction of origin of periodic surface error at 
surface I. 

One- dimensional array (101). Displacement in 
y direction of origin of origin of periodic surface 
error at surface I. 

One- dimensional array (101) for displacement Ax 
of coordinate system at surface I. 

One-dimensional array (101) displacement Ay of 
coordinate system at surface I. 

One- dimensional array (101), displacement Az of 
coordinate system at surface I. 

One- dimensional array (101) for angular rotation 
or tilt 6x in radians around the y axis at surface I. 

If THTAX(I) input greater than 2ir, it is assvimed 
to be degrees. 

One-dimensional array (101), angular rotation or 
tilt 6y in radians (or degrees if 0y > Ztt) around the 
X axis at surface I. 

One- dimensional array (101), angular rotation in 
radians (degrees if 6z > Ztt) around the z (optical) 
axis at surface I. 

One- dimensional array (101), semiaperture of 
obscuration in x direction at surface I. 

One- dimensional array (101), semiaperture of 
obscuration in y direction at surface I. If both 
BETAPX and BETAPY are entered, the obscuration 
is rectangular. If only one is entered, the obscuration 
is circular (i.e., BETAPX = 12., BETAPY = 0.,). 
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BETASX(I) 


VZER0(I) 


VPI(I) 

GZER0(I) 


One- dimensional array (101), semiaperture of 
vignetting in y direction at surface I. 

If BETASX and BETASY are both nonzero, the clear 
aperture is rectang\ilar. If only one is nonzero, the 
clear aperture is circular. 

NOTE: Normally, only those rays are traced which 
lie inside BETASX, BETASY and outside BETAPX, 
BETAPY. If anegative sign is attached to any of 
the four values at surface I, the logic is reversed 
so that only rays are traced which fall inside 
BETAPX, BETAPY or outside BETASX, BETASY. 
This produces the effect of an annular obscuration. 

One- dimensional array of length (2). Allows repre- 
sentation of percent vignetting of upper edge of 
entrance pupil. Expressed as fraction so that 
VZER0 = 0.5 reduces clear aperture of upper edge 
to BETAO/2. VZER0 = 0. or 1. represent no 
vignetting . 

VZER0(1) — Received signal (or transmitter) 

VZER0(2) L.O. 

Default: VZER0(1) = 0., 0., 

Same as VZER0 but for lower edge of entrance 
pupil. 

One- dimensional array (2). Entrance pupil 
obscuration parameter. Distance from center of 
entrance pupil to edge of upper obscuring aperture. 


33 



GPI(I) 

UZER0(I) 


VLAMDA 

NFRS 


GZER0(1) -* Received signal or transmitter 
GZER0(2) - L.0. 

Default: GZER0 (1) = 0., 0., 

One- dimensional array (2). Same as GZER0 but for 
lower edge of obscuring aperture. 

One- dimensional array (2). Radius of obscuring aperture. 
Arc of radius UZER0 passes through GZER0 (concave 
downward) and through GPI (concave upward). No rays 
are traced below the arc through GZER0 or above the arc 
through GPI. For a circular, centered obscuration 
UZER0 = GZER0, UZER0 must be nonzero if GZER0 
and/or GPI are entered. 

UZER0(1) — Received signal or transmitter 
UZER0(2) - L.O. 

Default: UZER0(1) = 0., 0., 

Spectral wavelength of operation of system being 
analyzed. All dimensional data (TS, BETAO, etc.) must 
be in same units as wavelength. 

Default: VLAMDA = 0.010611385 (mm) 

Integer input specifying number of intervals or output 
points for which spread function (ASF, PSF) arrays are 
computed. Dimensions of spread function array 
= (2- NFRS- 1) X (2- NFRS- I) • 

NFRS <51. 

Defavilt: NFRS = 26 

(giving spre^id function arrays = 51 x 51). 
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DELFS 


NXY 


PZER0(I) 


Spread function interval size. Nominally l/lO the Airy 
disc radius. Detector parameters are specified as inte- 
ger multiples of this value. 

Default: DELFS + °i^ : ‘BETllj" 

A 

where FL and BETAO are respectively the focal length 
and semiaperture for the received signal or transmitter 
optics. 

Integer input specifying number of grid points over 
entrance pupil semidiameter. The total pupil array 
is (2-NXY + 1) X (2-NXY + 1) 

NXY < 50 

Default: NXY = 25 

(giving pupil array = 51x51). 

Receiver — PZER0(1) = Radiant power density at entrance 
pupil of received signal optics. If PZER0(1) input 
negative, PZER0{1) = total radiant power incident 
(uniformly) on entrance pupil. PZER0(2) = total radiant 
power in trxincated gaussian beam of L.O. laser. 
Transmitter — PZER0(1) = total radiant power in 
truncated gaussian transmitter laser beam. 

Default: PZER0(1) = l.,l., 
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ZC^MEGA 

IPRi&Ga) 

IPR(&G(1) 

IPR(Z)G(2) 

IPR<&G(3) 




IPR<Z)G(4) 

IPR«5G(5) 

IPR«J)G(6) 


2 

Radius of 1/e intensity point of transmitter or L.O. 
gaussian laser beam. 

Default- transmitter: Z0MEGA = BETAO(l) 

Default- receiver: Z0MEGA = BETA0(2) 

One- dimensional array of integer flags for various 
program options. 

Default: IPR0G(1) = 1,0, 0,0, 0,0, 

= 0 - Paraxial analysis only 

= 1 - compute receiver or transmitter quality criteria 
= 0 - Pass 

= 1 - Output receiver ASF arrays or transmitter PSF 
array. 

= 0 - Compute receiver quality criteria with detector 

centered on optical axis, received signal and L.O. 
beams shifted due to effect of input beam obliquity, 
misalignments, IMG error, etc. 

= 1 - Center received signal and L.O. chief rays on 
detector 

= 2 - Center respective peaks of ASF for received 
signal and L.O. on detector. 

= 0 - Pass 

= 1 - Print un-normalized transmitter PSF 
= 0 - Pass 

= 1 - Print OPD (w(u, v)) arrays 
= 0 - Pass 

= 1 - Print pupil function 
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NDET 


IDETEC(I, J) 
IDETEC(I, J) 

IDETEC(1, J) 

IDETEC(2, J) 
IDETEC(3, J) 

IDETEC(4, J) 

IDETEC(5, J) 


Integer input specifying number of detectors to be 
analyzed. 

Defatilt: NDET = 1, 

Two-dimensional array (5, 5) defining parameters of 
detectors to be included in analysis. 

*— a } b I c I d J e } 

I = a,b,c,d, e, 

J = Detector number 

= a = circular, rectangular flag 

1 - circular detector 

2 - rectangular detector 

= b = Integer input specifying number of grid points 
across x dimension of detector. Total width Sx of detec- 
tor in X direction is Sx = 2* b* DELES 

= c = Integer input specifying number of grid points 
across y dimension of detector. Total width Sy of detec- 
tor in y direction is Sy = 2* c * DELES. Eor circular 
detector c = 0 and Sx = circular diameter. 

= d = Integer input specifying number of grid points 
by which detector center is to be displaced from the 
optical axis in the x direction. Detector displacement 
= Dx, and Dx = d* DELES. 

= e = Integer input specifying number of grid points by 
which detector center is to be displaced from the optical 
axis in the y direction. Detector displacement = Dy 
and Dy = e* DELES 
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Default: 


a = 1 
b = 10 
c = 0 
d = 0 
e = 0 

NOTE: Any combination of b and d; c and e; or off-axis 
chief ray shift causing the detector to fall outside the 
spread function array (NFRS) will abort the run. 

FLAG An integer flag used to reinitialize data. Primarily 

used between consecutive runs to wipe out all data from 
a previous run. This is submitted as an independent data 
deck with title card, ^INP card with FLAG = 1, and 
S^END card. Stacking this between two data decks 
accomplishes the initialization for the second set of data. 

2.7 INTERPRETATION OF OUTPUT DATA 

The amount of output data printed is controlled by the IPR0G 
flags. The system parameters as input will always be printed; the 
paraxial ray traces and computed first order parameters will always be 
printed. For the full transmitter or receiver analysis the full set of 
quality criteria will always be printed. 

a. Paraxial Data 

• Entrance pupil position T(0) = , the value printed here 
is the distance from the paraxial entrance pupil to 

the first surface. A negative value indicates that the 
- entrance pupil lies to the right of the first surface. 

• TEXIT — Paraxial exit pupil position, TEXIT = distance 
from surface (N-1) to paraxial exit pupil. 

• XLINV — Inverse object distance as input or as computed 
from input of MAG. 

• FOCAL LENGTH — computed paraxial focal length. 
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BACK FOCAL LENGTH — computed paraxial back 
focal length = distance from surface (N-1) to paraxial 
focal plane. 

TS (N- 1) — Spacing actually used in analysis as 
distance from surface (N- 1) to image or detector 
surface. Depends on IFLGl, TS(N), etc, 

FL- when XLINV 4 0. FL = MAG/XLINV so that 

an input wavefront radius at the entrance pupil of 
R = L = 1 /XLINV will give an output wavefront 
radius <?f FL. 

Chief Ray Data 


Field Angle — is the value input as H 

Chief Ray Coordinates — x, y, z are the coordinates 
of the intersection of the chief ray with the image 
surface. 

Direction Cosines — 1, m, n for chief ray in image space 

RW — Radius of reference wavefront to be used to 
compute OPD values over exit pupil. 

TR — Location of reference wavefront, TR = Axial 
distance from surface (N-1) to reference wavefront. 

For axial case TR = TEXIT. 

XTILT — For off-axis case, reference wavefront 
may be tilted or rotated XTILT = THTAX for refer- 
ence wavefront. 

YTILT — = THTAY for reference wavefront. 

XDISP— with tilts and decentrations, the chief ray 
may become a skew ray so that it does not intersect 
the optical axis. XDISP is the distance from the chief 
ray to the optical axis at its'point of closest approach. 

XDISP = DLTAX for reference wavefront. 

YDISP— Same as XDISP, YDISP represents distance 
from optical axis to chief ray at point of closest 
approach. 

YDISP = DLTAY for reference wavefront. 



c. Wavefront Error Information 

After all the OPD (W(u, v)) values have been computed 
over the reference wavefront, wavefront statistics are computed. 

SI 

u V 


=2^W{u.v) 

U V 



• Maximum un-normalized error = WMAX 

• Minimum un-normalized error = WMIN 

• Maximum normalized error = WMAX/K 

• Minimum normalized error = WMIN/x. 

WMAX = maximvim value of W(u, v) 

WMIN = minimum value of W(u, v) 

• Approximate Strehl ratio = I(cr) 

I(o-) = 1. - (2 -t:-STD DEV)^ 

d. Spread Functions — (ASF and PSF) IPROG(2) = 1. — 

For the receiver analysis, selection of this option prints out the arrays 
FR and FI for the received signal and the local oscillator (ASF(x, y) = 
FR(x,y) + iFI(x,y)). 

labelled 

"ASF (REAL), FOR RECEIVED SIGNAL” (FRl) 

"ASF (IMAGINARY), FOR RECEIVED SIGNAL” (FIl) 
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"ASF (REAL), FOR LOCAL OSCILLATOR" (FR2) 

"ASF (IMAGINARY), FOR LOCAL OSCILLATOR" (F12) 

"PHASE" is also printed where PHASE = ARCTAN (FI/FR). 

The ASF values are normalized to 100 for compactness and 
the normalizing factor is printed. The ASF arrays will be (2-NFRS-l) 
X (2*NFRS-1) up to 72 x 72 beyond which only the central 72 x 72 points 
will be printed.'- If the array is smaller than 36 x 36, the entire array 
will be on a single page. Arrays larger than 36 x 36 will be printed 
one quadrant per page. The quadrants are related as per Fig. 2-5 
with the 0, 0 point in the upper left hand corner of quadrant 4. 


3079-4 


2 

1 
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V 

0,0 
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Fig. 2-5. 

Output quadrants for ASF, PSF, pupil 
function. 


PSF — The PSF or intensity value for a point in an array is 
PSF(x, y) = (FR(x, y))^ + (FI(x,y))^. The PSF array is printed for 
the transmitter. All the comments for the ASF printout apply here. 

INTERVAL — This is the spacing between adjacent ASF or PSF 
points. INTERVAL = DELFS. 

NUMBER INTERVALS = 2'NFRS-l 

e. IPR0G(4) =1. — This option causes the printing of the 
floating point tabulation of the transmitter PSF (intensity) values. 

f. IPR0G(5) = 1. — This option causes the floating point 
tabvilation of the OPD values to be printed. 
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g. IPROG(6) = 1. — This option prints the complex 
components of the pupil function normalized to 100. As with the ASF 
and PSF, the output may include only the central 72 x 72 of a larger 
array. 

2.8 COORDINATE - SIGN CONVENTION SUMMARY 


The coordinate system for ray tracing is a local coordinate sys- 
tem where the vertex of the surface being traced is the origin. The 
spacing TS combines with any misalignment parameters to transfer 
the coordinate system from surface to surface with the ray. Some 
key points regarding sign conventions and coordinates are summarized 
here. 

• Surface parameters — A surface parameter with a 
positive sign causes the surface to be deviated to 
the right of the tangent plane, hence the negative 
sign indicates a deviation to the left of the tangent 
plane. Some of those parameters are: 

RH0S, RH01, RH02, RADS, DF0RM, AMP. 

• Spacing — TS. A positive value for TS(1) specifies 
that surface I + 1 lies to the right of surface I. 

• Finite object — nonuniphase input beam. A positive 
value for X LI NV indicates that the entrance pupil 
lies to the right of the object point. The input wave- 
front is thus divergent. For MAG input positive for 
a system with a positive focal length, the input 
wavefront is divergent and the exiting wavefront is 
convergent. 

• H-input beam obliquity — A positive H indicates an 
input beam incident on the entrance pupil from below 
the optical axis. 

• Misalignments - DLTAX, DLTAY, DLTAZ, THTAX, 
THTAY, THTAZ- refer to Section 4. 6. 

• Detector — The optical axis is the reference here. 

For a centered system with no input beam obliquity, 
the detector, the received signal and the L. O. beam 
will be centered on the optical axis. If there is an 
input beam obliquity or misalignment in either the 
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received signal or the local oscillator, that beam 
will be shifted from the center of the detector 
while the other beam remains centered. The 
detector will always remain centered regardless 
of the shift in the beams unless the IPR0G(3) 
option is invoked or the detector is displaced 
via IDETEC(4, J) and/or IDETEC{5, J). Notice 
that IPR0G(3) centers the spread functions on 
the detector while IDETEC(4 and 5, J) shift the 
detector with respect to the optical axis without 
changing the location of the beams. 


2.9 SAMPLE CASES 


Computed results are included for four sample cases. There 
are two transmitter cases and two receiver cases. The descriptions 
of these cases are as follows: 

a. ’'TRANSMITTER TEST CASE", Appendix A - This very 
simple case is an example of an analysis using a "perfect surface". 
Surface No, 1 is a dummy surface, surface No. 3 is the image surface 
and surface No. 2 is the perfect imaging surface. For this perfect 
system, the focal length and back focal length are equal at 100 mm. 

The transmitter quality criteria are determined at the (Fraunhofer) 
focal plane. The overall transmitter efficiency is 92% rather than 
approaching 100% because of the gaussian beam. 

An examination of the data deck for this case would indicate 
the presence of data for a complete receiver. The transmitter analysis 
results from setting N(2) = 0. 

b. "TRANSMITTER TEST CASE No. 1", Appendix B - This 
system consists of a meniscus lens represented by surfaces No. 2 and 
No. 3 and a paraboloidal mirror represented by surface No. 5. Sur- 
face No, 1 is a dummy located at the position of the laser input; surface 
No. 4 is a dummy at the mutual focus of the meniscus and paraboloid. 
The combination of the meniscus and paraboloid produces a perfectly 
collimated beam so that surface No. 6 is a perfect surface chosen to 
focus the beam at a distance of 8 km (back focal length — 

8.00178987* 10^ mm). The laser beam defaults to a total power of 
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2 

1 W. The l/e radius is specified to be 0. 925 mm (Z0MEGA = 0.925) 
and the far field intensity is sampled in increments of 50 mm (DELFS 
= 50. 0). The reference wavefront has a radius of 7. 9 km 
(RW = 7. 93148013*10^ mm) and is located 70 m to the right of surface 

4 

No. 6 (TR = 7. 03097397* 10 ). The quality criteria are essentially 

self-explanatory with the added comment that the apparently low value 

for overall transmitter efficiency is due to the disparity between the 
2 

l/e radius (0.925 mm) and the truncation radius (BETAO = 2.25). 

c. "RECEIVER TEST CASE No. 1", Appendix C - The 
received signal optics for this system are for a real system. The 
local oscillator optics is a perfect surface of aperture and focal length 
to give an Airy disc about five times the diameter of that for the 
received signal. The peak value for the ASF(REAL) for the received 
signal (FRl) is 1304. 5 while the peak value for the ASF (Imaginary) for 
the received signal (FIl) is 600. 1 indicating that the peak intensity is 
2, 062. 10^ and that the focus can be improved (to reduce FIl). The 
ASF values for the local oscillator indicate good focus so that the peak 
intensity will be 2. 004^ = 4. 0160. 

The ASF arrays are output for the two beams allowing an examina- 
tion of the respective distributions. The FI arrays are converted to 
phase and exhibited since default conditions are relied on for the quality 
criteria. They are based on a circular detector of diameter equal to 
that of the Airy disc for a system of clear, uniformly illuminated 
aperture, and relative aperture = 8. 0. 

d. "RECEIVER TEST CASE No. 2", Appendix D - This data 
is for a typical set of received signal and local oscillator optics. The 
quality criteria are determined for the default case plus three other 
detectors: 

(1) Detector No. 1 — Default case — circular, centered 
detector with diameter equal to the Airy disc. 

(2) Detector No. 2 — Square, centered detector with 
sides equal'to the diameter of the Airy disc. 
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Detector No. 3 — Circular detector with same 
dimensions as Detector No. 1 but shifted from the 
beam center by a distance equal to the radius of 
the Airy disc. 

Detector No. 4— Square detector with the same 
dimensions as Detector No. 2, but shifted from 
beam center by the Airy disc radius. 



SECTION III 


SAMPLE SYSTEM ANALYSIS 


The system chosen for exercising this computer program is the 
optical mechanical subsystem (OMSS) developed for NASA xinder con- 
tract NAS 5-21859. This package is intended as an engineering model 
of a spaceborne data link and is one of the most advanced systems 
designed specifically for heterodyne communications work. The system 
has been designed and analyzed using existing optical computer programs 
and represents a nearly perfect optical system. As will be seen, the 
excellent performance predicted by a conventional program, ACCOS V 
is confirmed by LACOMA. Since the performance is nearly perfect, 
the differences between LACOMA and conventional programs are masked. 
An example of a less perfect system is presented at the end of this 
section to illustrate the advantages of LACOMA over programs that 
optimize intensity. 

The OMSS consists of two focusing reflective surfaces, a para- 
bolic primary and an elliptic secondary. This particular configuration 
provides essentially perfect performance when properly aligned, and, 
even in the presence of reasonable tilts and decentrations, performance 
is still excellent. 

The analysis was carried out in three steps. First, a reference 
case using the same size optics with the same obscurations but with a 
"perfect" surface in the receive and the L. O. paths was analyzed to 
provide a measure of the best possible performance vinder the given 
restraints. The configuration of this system is shown in Fig. 3-1; the 
optical input data and the results are listed in Table 3-1. 


Preceding page blank 
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TABLE 3-1. 


Reference System Specifications and Results 


Specifications 




Receiver F/No. 

F/8 



Entrance Pupil Diam 

16.51 cm 



Obscuration Diam 

3.20 cm 



L. O. F/No. 

F/40 



Detector Diam 

0. 199 mm 



LACOMA Results 




Phase Match Efficiency 

- 

1. 000 

-0.00 dB 

Optical Transmission 


0.964 

-0. 16 dB 

Direct Detection Efficiency 


0.743 

-1.29 dB 

L.O. Illtimination Efficiency 


0.823 

-0. 84 dB 

Theoretical Antenna Gain 



93. 78 dB 

Receiver Efficiency to I. F. 


0. 504 

-2.98 dB 

Effective Antenna Gain 



90.80 dB 


T1200 


The second step of the analysis was to evaluate the system 
without alignment errors to compare it with the reference case. Four 
configurations were evaluated; an on-axis case with and without a beam- 
splitter that will be used to add a beacon to the system, and an off- 
axis case also with and without the beamsplitter. The off-axis angle 
was chosen to equal the maximum angle encountered during the signal 
acquisition phase. The third step of the analysis was an evaluation of 
system sensitivity to manufacturing tolerances. 

The optical configuration for the second and third parts of the 
analysis is shown schematically in Fig. 3-2. Several dummy surfaces 
were included to allow restoration of axes and to provide means to 
perturb optical elements without changing the basic constants of the 
system. The surfaces and their functions are listed below: 

' (1) Parabolic primary mirror 

(2) D\immy used to restore tilts and displacements of 

the primary mirror, or to add spacing errors between 
the primary and secondary. 
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RECEIVER SCHEMATIC 




Fig. 3-2. OMSS optical schematic. 
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(3) Elliptical secondary mirror 

(4) IMG mirror 

(5) Dummy to restore the optical axis to 90° after 

reflection from the IMG. The sum of the tilt 
angles of surfaces 4 and 5 must always equal 
90°. 

(6) Aperture stop 

(7) First beamsplitter surface 

(8) Second beamsplitter surface 

(9) Dummy to restore tilted axis introduced by surface 
No. 7 

(10) Dummy to restore offset between optical axis and 
chief ray caused by the beamsplitter 

(11) Detector plane for receiver path 

(12) L. O. perfect surface and aperture stop 

(13) Detector plane for L.O. 

The constants to describe each surface for the nominal on-axis 
case are listed in Table 3-2. These values were determined during the 
OMSS design phase, and have been optimized for best Strehl ratio 
(diffraction focus). Note that there are three reflective surfaces in 
the system; both the refractive index and the direction of propagation 
(spacing) reverse sign at these surfaces. 

Local Oscillator — At the time of this analysis, there was no 
final design data for the L.O. Accordingly, to avoid introducing 
degradation from this source, the basic design data was used to gen- 
erate a "perfect" L.O. The planned L.O. will operate as a f/40 system 

2 

with a gaussian beam distribution truncated at the l/e points, and 
with no central obscuration. A perfect L.O. that meets these criteria 
is also shown in Fig. 3-2, and is specified in Table 3-2. 

Figure 3-3 shows a printout of the data deck as it appeared after 
being input to the computer. Additional parameters shown in Fig. 3-3 
are defined in Section II. 
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0’-1 SUBSVeTL'i, .'JO OFFSETS OR TILTS 
JIMP 

PETA0(1)=^^*-'.^3»5. , f’ZEROd) =-i. ,1., I Ft G 1 ( 1 ) = 1 , IPROG ( 1) = 1 , 0 , 2 , 

N( 1) =-ll,-2i X*1iJS(ll = -1. , -1. ♦l.,3»-l.,-4.03 06 20,-l.,3*-l.,M(H=6»12 
SKAPA (1) = u . » 1 . , . 7173F , o i T A^ X ( l ) =15 . , 

R ADS ( 1) = “^lc J . , C . , 141 . d :i7 3, 5*'!. , 

TS (1) = C 2. 9 J98F , 13 3. 4 + 557, i . ,-12. 425,-2 . , u . , -54. 5353 , 2*0 , 

TS(121=4Cu. ,l . , RAC‘S( 121 . , J . , 3RS(12) = -1., XMUS ( 12 ) = 2'^1 . , 

DELFS = 4. l‘+1657E-3, NQET = 1, I OET EC ( 1 , 1 ) = 1 , 2 4 , 3 ♦ u , 

THTAY (4) =2'^- . 7 55 39ftlT , 

THTAY(7)=.73539315,0.,-.73539515, 9LTAY(1:)=.25396, 

?£NO 


Fig. 3-3. Input deck for OM system. 
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The results for the runs without alignment errors are shown in 
Table 3-3. The results show that the aberrations caused by off-axis 
operation and the astigmatism caused by the beamsplitter are both 
negligible. (In this table the overall efficiency to the i.f. is shown in 
dB degradation below the theoretical maximum antenna gain, and is 
also normalized to the resvilts of the reference case. ) 

For the third step of this analysis, the distortions introduced 
into the OMSS were selected to allow comparison with the analysis 
conducted using ACCOS V as listed in the design report for the OMSS. 
Four perturbations were analyzed: 

(1) The primary- secondary spacing was varied 

(2) The detector was displaced from the focal plane 

(3) The primary mirror was transversely displaced 

(4) The primary mirror was tilted. 

The results are presented graphically in Figs. 3-4 through 3-7. 
For the first two cases where the evaluation is conducted on- axis, the 
results are directly comparable to those obtained using ACCOS V. The 
plots show Strehl ratio, normalized heterodyne efficiency, direct 
detection efficiency, and local oscillator efficiency. There are several 
interesting results: 

(1) Because the local oscillator has been spread out by 
the f/40 optical system, the amplitude profiles of 
the L.O. and received signal tend to become more 
closely matched as the received signal is defocused, 
thereby increasing the local oscillator efficiency 
and partially compensating the loss in signal power. 

(2) Heterodyne efficiency is not as sensitive to focus as 
is the Strehl ratio, partially because of the above 
observation. 

(3) The region around the focal plane of the detector 
would be symmetrical in the absence of aberrations. 

The aberrations that are present cause a slight 
asymmetry in the Strehl ratio, which is barely 
evident on the curves. However, the heterodyne 
efficiency shows amuch greater asymmetry, 
implying that small aberrations affect the phase 

of the signal more than the amplitude. 
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Reference System 
OMSS-No Beamsplitter 
On- Axis 
0.13° Off-Axis 
OMSS w /beamsplitter 
On- Axis 
0.13° Off-Axis 


TABLE 3-3. 

Results for Nominal OMSS Configuration 


Phase 

Match 

Efficiency 

Direct 

Detection 

Efficiency 

L. O. 

Receiver 

Efficiency to I. F. 

Efficiency 

Decimal 

dB 

Relative 

1. 000 

0. 743 

0.823 

0. 504 

-2,98 

1. 000 

0.996 

0.742 

0. 820 

0.496 

-3. 04 

0.985 

0.993 

0.735 

0. 820 

0.488 

-3. 12 

0, 967 

0.990 

0.741 

0.821 

0.495 

-3. 05 

0.983 

0. 990 

0.736 

0.824 

0.489 

-3.11 

0. 970 


T1202 
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Fig. 3-4. Primary- secondary separation error. 
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DETECTOR PLANE DEFOCUS, mm 


Fig. 3-5. Effects of detector position 
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PRIMARY DECENTRATION, mm 


Fig. 3-7. Effects of primary offset. 
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For cases 3 and 4 where the primary has been tilted or 
displaced, the results cannot be compared directly with previous results 
due to differences in methods used to evaluate the system. For the 
L.ACOMA r\ins, the input beam was introduced along the undisturbed 
optical axis and the IMC was adjusted to place the chief ray of the 
received signal onto the detector center as would occur during the track 
mode in the operating system. The values for IMC correction were 
determined by paraxial ray trace using a desk calculator. The LACOMA 
option to place the peak intensities of the two beams on the center of 
the detector was selected to eliminate residual errors in IMC position 
and to compensate for peak shift due to off-axis aberrations. An alter- 
native approach leaves the IMC centered and compensates for the 
angular displacement by evaluating a received signal from an off-axis 
location. The aberrations caused by the first approach are more 
severe, as is evident if the Strehl ratios in Figs. 3-6 and 3-7 are com- 
pared with corresponding figures in the OMSS design report. 

For these cases, as shown in Figs. 3-6 and 3-7, the heterodyne 
efficiency falls off more rapidly than the Strehl ratio. This result is 
not unexpected, as the heterodyne signal is more sensitive to wavefront 
tilt, thus the phase match efficiency suffers. Also of interest is the 
L.O. illumination efficiency which decreases for these cases, thus 
adding to the degradation. 

Overall, the results of this analysis tend to confirm the con- 
clusions reached in the OMSS design report. The heterodyne signal 
has been shown to be more sensitive to alignments that cause angular 
errors than those that cause only defocusing. This result is beneficial 
to the OMSS, since the primary- secondary spacing is the most difficult 
parameter to control. 

Other Analyses 

The LACOMA program has been used to analyze several other 
problems of interest; two of these are summarized here. 
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Field of View — The nominal field of view of a heterodyne 
receiver has been generally considered to be k\/d, where d is the 
receiver aperture diameter, and k is a constant dependent upon the 
illumination factor, detector size and definition of field of view. To 
evaluate k for a round detector equal in diameter to the Airy disc, 
LACOMA was used with a perfect system and a series of detectors 
spaced sequentially at varying distances from the optical axis.. 

Figure 3-8 shows the results for a uniform L.O. illumination. From 
this curve, the full field of view to the 3 dB points implies that k should 
be 1.38 for the uniform case. 

Germanium Galilean Telescope — A four-element telescope 
configuration meant for use in a lO.b-pm heterodyne system is shown in 
Figure 3-9. This telescope was designed for maximtim Strehl ratio, 
using existing design programs. {This particular telescope has suf- 
ficient spherical aberration in this configuration to reduce the Strehl 
ratio to about 0. 7. ) Since this system was designed as an afocal sys- 
tem, a perfect surface was introduced into the afocal beam to provide 
a focal plane for evaluation, and a perfect f/40 L.O, was also selected. 
The lens system was evaluated, using the predicted optimum spacings, 
to determine the degradation caused by the spherical aberrations; also, 
with the spacings between the first and second elements increased by 
0.25 and 0.5 mm to ascertain the effects of spacing tolerances. The 
results, shown in Figure 3-10, were unexpected; the overall efficiency 
increased from -4. 55 dB to -2, 58 dB, which is better than the -2. 74 dB 
expected from a perfect system with a uniform L.O, This result can be 
explained by reference to the phase match efficiency and the L.O. 
illumination efficiency. Both of these values have increased sig- 
nificantly from the nominal position, even though the total signal power 
has dropped. This means that at a location significantly away from the 
best energy focus, the wavefront is nearly plane and the energy dis- 
tribution is more \iniform. Thus, in the presence of significant aber- 
rations, optimum performance for a heterodyne system does not neces- 
sarily correspond to that for a direct detection system. LACOMA is 
the first program to put a quantitative value on this phenomena. 
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Fig. 3-8. Field o£ view variation 
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Fig. 3-9. Optical schematic for afocal evaluation. 
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Fig. 3-10. Results on Galilean telescope. 
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SECTION IV 


DETAILED PROGRAM DESCRIPTION 


Upon reading the input data, the program determines whether 
it is to execute a transmitter analysis or receiver analysis based on the 
input of one or two sets of optical data. Specifically, it checks whether 
values are entered for N(l) and N(2) which identify the number of sur- 
faces in the two sets of data for received signal optics and local oscil- 
lator optics, respectively. If N(2) = 0, then a transmitter analysis is 
performed. In addition to the surface number parameter, many of the 
other system parameters are two-element arrays; e.g., M(l) = surface 
number of aperture stop for received signal (or transmitter) optical 
train while M(2) = surface nvimber of aperture stop for local oscil- 
lator optical train as is the case for H(l), H(2), FL(1), FL(2), etc. 

As most of the computations to be discussed are performed identically 
for transmitter, receiver, and local oscillator, the subscript notation 
will be omitted except where confusion might arise and the distinction 
between transmitter and receiver will be omitted except for those 
computations unique to a given case. 

4. 1 FIRST ORDER PARAMETERS 

Two paraxial rays are traced through the optical trains to deter- 
mine the first order parameters of the particular optical train. These 
parameters are: 

a. Paraxial entrance pupil position T is the distance 
from the paraxial entrance pupil to the first optical 
surface. 

b. Paraxial exit pupil position TEXIT or TR is the 
distance from the last optical surface (surface N- 1) 
to the paraxial exit pupil. 
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d. 


Inverse object distance XULNV or h is the inverse 
of the distance from the object plane to the paraxial 
entrance pupil. 

Focal length 


e. Back focal length, which is the distance from the (N- 1) 
surface to the paraxial focal surface. 

f. TS(N-l) is based on the input data and control param- 
eters, and is the spacing actually used in the analysis 
as the back focus. 

g. Paraxial Computations are where the two paraxial 
rays are traced with starting values at the first sur- 
face of the optical train. 

pQ = input 

b = 0. 
o 

o o 

a = 0 . 

o 

With this data as input the two paraxial rays are traced using 




b = b 
s s- 


s s- 


+ b , t , / p 
s-1 s-F '^s-1 

(4-1) 

^ ^^s- 1 Ps Ps 

(4-2) 

+ a 1 T , / M , 

s-1 s-1 ^^s- 1 

(4-3) 

^^"s-l ■ ^"s^ “"sPs 

(4-4) 


The parameters of these paraxial rays are illustrated in Fig. 4-1. 
The data for the paraxial rays at the first surface, the aperture stop, 
and at the rear surface are then used to compute the desired first 
order parameters. This data is designated as follows: 
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Fig. 4-1. General paraxial ray notation. 
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(B , b , O' , a ) at the first surface 
o o o 

(p^. b » Q , a ) at the aperture stop 

m m m ^ 

(p I b , O' , a ) at the rear optical surface 
“r r r r 

(The r subscript denotes the next to last or (N- l)st surface) 



( 4 - 5 ) 


(o'™ P - O' p ) 
m ^^r r ur 


EXIT 


^r P: 


- O' b 
m m r 


( 4 - 6 ) 


The focal length F and back focal length BF for infinite con- 
jugates are computed: 


, F 



( 4 - 7 ) 


BF 


Pr^r 

b 


( 4 - 8 ) 


If a finite conjugate distance is entered in the form of 
(XLINV) to introduce a quadratic phase component into the input beam, 
the focal length and back focal length are recomputed; 


F 


b 

r 



+ a 


r 




BF 



+ O' L 
r o 



b 

r 


+ a 

r 




( 4 - 9 ) 


( 4 - 10 ) 
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Optionally, a value may be specified for magnification as the 
ratio of the input beam divergence versus the output beam convergence: 

-b 

Magnification = mag = - (4-11) 

r 

from which the inverse object distance is computed for use in the other 
calculations, 



o 


-mag b^ 

(1 + mag p^a^) 


(4-12) 


Notice that for the finite conjugate case, BF is the location of the 
paraxial image for the given object position and the focal length FL is 
not the classical paraxial focal length but is 


FL = mag • LO, 


(4-13) 


This value is used in the computation of the input ray angles 
for tracing rays. 

4.2 IMAGE LOCATION 


The actual value used to locate the image surface of intensity 
plane is subject to several options. The program computes the parax- 
ial back focus BF based on the specified object location. There is a 
value input as tj^ ^ for the spacing between the last ((N-l)st) surface 
and the image (detector or intensity) plane. The option is available 
for specification of which of these is to be used in the analysis. The 
parameter IFLGl controls this option. 

^k-1 =BF if IFLGl =0 (4.14, 
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tj^_ 1 = tj^_ 2 if IFLGl 1 (4-15) 


where 

I 

2 is ih® value to be used in the analysis. 

Additionally, the value tj^ which is usually ignored since it has 
no meaning, can be used to specify an incremental shift of the image 
location. Thus, the value of tj^ is always added to the back focus or 
image distance. 

tj ^_2 = iFLGl =0 (4- 14a) 

^N-1 = ^N-1 = i (4- 15a) 

For example, if it is known that the best focus for a given sys- 
tem is 0.001 from the paraxial focus, the program can compute the 
paraxial back focus very accurately and then shift by 0. 001 to the posi- 

I 

tion of best focus. The value of 2 thus determined is printed out as 
TS (N-1). 

4.3 RAYTRACE-OPD 

The LACOMA ray trace constitutes an optimal combination of 
accuracy and speed of computation. The ray trace equations are the 
result of many years of development, and permit the trace of both 
meridional and skew rays through any surface which is continuous and 
analytic . 

Notation — (Figure 4-2). A given ray is described by the 

parameters X , Y , Z , 1 , m and n . X , Y and Z are the 
^ sssss sss s 

coordinates of the point of intersection of the ray with surface ”s". The 
origin of the coordinate system is at the vertex of the surface. The 
Z axis is the optical axis, the Y axis the meridional axis and X the 
skew axis. The optical direction cosines of the ray exiting surface ”s" 


70 




are 1 , m and n . In a medium following surface "s”, the refractive 
s s s 

index is u and 
s 


i2, 2.2 2 

1 + m + n = u (4- 16) 

s s s s 

Surface "s" itself is specified by some sagitta function ^ = ^(Q) and, 
for a general rotationally symmetrical asphere, 

P Q ^ 

„ so s 

® 1 + (l- Kp ^ Q ^ 

\ '^so s 

(4-17) 


4 , 7T 


6 . - 


T72 


+ »Q +pQ +YQ 

• o ey 


+ 



+ Y 


2 

s 


p = vertex curvature of surface = 1/R 
so so 


etc. = general aspheric deformation coefficients 
(not related to o- or p for paraxial ray 
trace) 


K 


s 


(l-€ 



conic coefficient of surface 


Conic sections — If the conic section surface is thought of as a 
mirror with two aplanatic foci and F 2 as in Figure 4-3, the conic 
eccentricity can be related to the distances Sj and S^, respectively, 
measured from the first focus F^ to the mirror and from the mirror to 
F^. In Fig. 4-3, is positive and is negative. The vertex curva- 
ture for the conic section is 




(4-18) 
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Fig. 4-3. Conic eccentricity parameters. 
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and the eccentricity is 


e = 


^^2 
^ 1 - ^2 


(4-19) 


Eccentricities for some typical surfaces are 


Surface 

Sphere 

Paraboloid 

Ellipse 

Hyperboloid 


0 . 

1 . 

0 < € < 1, 

« > 1 . 


K 

1 . 

0 . 

0 < »< < 1 
K < 0. 


The sagitta expression for a sphere would be 


^ = 


P Q 
^so s 


1 + (l - P ^ Q 

\ so s / 2 


(4- 17a) 


and for a paraboloid 


; = 


P Q 
so s 


(4- 17b) 


a. Ray Trace Equations — Denoting surface "s" with the 

"s" subscript and the preceding surface by the "s-1" subscript. 
Surface-to-surface transfer equations are 


X = X , + (t , - Z 1 + Z ) 

s s-1 n , \ s- 1 s-1 s/ 

s- 1 


(4-20) 


Y = Y , + — ^ (t , - Z , + Z ) 
s s-1 n 1 \ s- 1 s-1 s/ 

s-1 


(4-21) 


74 



Refraction at surface ”s" converts optical direction cosines 

1 1 , m , , and n , into 1 , m , n . 

s-1 s- 1 s-1 s s s 


2 = f , + P X (n , - n ) 

s s-1 s ss-1 s 


(4-22) 


m = m ,+P Y (n ,-n) 
s s-1 s's-1 s' 


(4-23) 


or 


i = H - p X n 
s so ^s s s 


(4-22a) 


m = m - p Y n 
s so s s s 


(4-23a) 


where 


i = i , + P X n , 
so s-1 s s s-1 


m = m + p Y n , 
so s-1 s ss-1 


(4-24) 

(4-25) 


and 


p = the inverse of the subnormal 
s 


I iL 

Q dQ 


so 


1 -7k P 2 q2 U 

Is so s k 


4.S Q, 4 6 4 8Y^ Q;' 4 


(4-26) 


For the general asphere the appropriate values of Z and n 

s s 

are found by iteration. Dropping the "s” subscript temporarily and 
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substituting i and i + 1 to indicate successive iterated values, the 
iterations on Z and n are 


1 - 


Z. p." (£ , X. + m , Y.) 

11 s-1 1 s-1 1 


n 


s-1 


1 + 1 


P." (e 1 X. + m , Y.) 
1 ' s-1 1 s-1 i' 


1 - 


n 


s- 1 


(4-27) 


2,2 .2 

p + n. -2. 


n 


m. 

1 


2 p n. (2. X + m. Y ) 

^s 1 ' 1 s 1 s' 


i + 1 


2 fn. - p " (X 2. + Y m.)l 

1-1 S SI S iJ 


(4-28) 


p. - 2 “ - m + n.^ (l + p ’• Q 

SO so 1 \ ^s s / 


zfn. fl + p"^Q^)-p"(£ X +m Y) 

1 1 \ s s / '^s ' so s so s J 


(4-29) 


using starting values of Z^ = 0 and n. = 


b. Spherical Surfaces — Though the above algorithms include 
spherical surfaces as a special case, the tracing through spherical 
surfaces is handled separately to expedite computing time. The trans- 
fer equations are 


X 


s 


X 

so 


+ H 

s 


^ 1 
s-1 


Y 

s 


Y 

so 


+ H m , 
s s-1 


(4-30) 

(4-31) 


i 
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Z = H n , 
s s s- 1 


'"so "‘s-l ‘ n , (^s-1 ^s-l) 

s- 1 


X = X , + 


Y = Y , + - (t j - Z ) 
so s-1 n ,\s-l s-1/ 

s- 1 


The refraction equations are 


n = n - - P P X 
s s-1 s s 


m = m , - P p Y 
s s-1 ss s 


n =n ,-p(p Z -l) 
s s-1 s\s s / 


where H and P are determined as follows: 
s s 

! 2 2 

(A) = p ( X + Y 
' ^s \ so so 


) 


(B) = n , - p (x £ 1 + Y m , ) 

' s-1 ^s \ so s-1 so s-1/ 


^^s-1 S ■ '^s-1 


2 

L^^s-1 


- P3 (A) 


H = 


(A) 


s (B) + p cos i 

3“ i S 


( 4 - 32 ) 

( 4 - 33 ) 

( 4 - 34 ) 

( 4 - 35 ) 

( 4 - 36 ) 

( 4 - 37 ) 

( 4 - 38 ) 

( 4 - 39 ) 

( 4 - 40 ) 

( 4 - 41 ) 
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11 cos r =11 
'^s s 

Ps = Pg COS cos (4-43) 

These algorithms are listed here in the form they are used in the com- 
putation which is necessary to preserve proper sign conditions. 

c. Special Surface Forms — Some surface forms which can 
be accommodated are noted below; most of them can be combined. 

Cylinders 

Toroids 

Normal surfaces with cylindrical error 
Periodic surface forms 

These surface forms, mentioned above, are described by the 
svimmation 


Z 2 . 

L , cos 1 
s- 1 s 


+ 1 


(4-42) 


S 

is the expression for a cylindrical surface of curvature whose 
axis may be rotated through an angle (p. Fig. 4-4. can combine 
with to generate a toroidal surface where the surface (curve) gen- 
erated by is rotated such that its vertex describes an arc of radius 
l/p^. The arcs of pj and p 2 lie in orthogonal planes. In Fig. 4-5(a) 
it can be seen that the cylinder is generated by p^ when (f = p^ = 0. 
Figure 4- 5(b) is the cylinder generated by p^ when p^ = ^ = 0. The 
surface (curve or arc) Cj may be the cross section of a sphere, any 
conic or higher order surface form. 

The periodic surface errors may be a radially symmetrical 
sinusoidal (cosine in this case) variation of sagitta, a radial sine func- 
tion or a combination of the two. These radial periodic errors may be 
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3079-17 



Fig. 4-4. Cylindrical surface due to RHOl but 
with ^ ss 0 . 


NAS 5-21898 

Final Technical Report 

Fig. 4-4 
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Fig. 4- 5(a). Cylindrical surface due to RHOl only. 
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Fig. 4-5(a). 
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j 



REPRESENTS: 

RH01 = 0 
RH02 * 0 
PHI =0 


;e due to RHOZ only. 


81 


displaced with respect to the axis by U and V in the X and Y directions, 
respectively. Figs. 4-6, 4-7, 4-8, 4-9, and 4-10. 

In the fornaulation of these surface sagitta, the following defi- 
nitions are used: 


Pj^ = curvature of cylinder (=1 /Rj) (RHOl) 

Py = toric rotation curvature or secondary cylinder 
curvature ( = 1/R2) (RH02) 

a = constant of rotational symmetry (CRS) 


b = cylindrical constant (CC) 


<p = angular orientation of 
cylinder (PHI) 


= angle between cylinder axis 
and X axis for 

= angle between cylinder axis 
and Y axis for P^ 


K = conic coefficient ( P^ surface) 


0, P, y, 
etc 


higher order (aspheric) surface coefficients (Pj^ surface) 


^1*^2’^ 3^ periodic surface error coefficients (AMP, FREQ, C3) 

U and V = displacement in X and Y of axis of radial symmetry of 
periodic surface error (USI, VSI) 


= xj + Y^ (4-45) 

2 2 2 2 2 2 
Q, = X (a + b sin t^) + Y (a + b cos o) + 2b X Y sincjcosus 

Is ' s s s 

(4-46) 
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Fig. 4-6 
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Fig. 4-7. Sinusoidal surface error with point synimetry of error decentered with respect 
to optical axis . 
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Fig. 4-8. Sinuso 
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Fig. 4-9 
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Fig. 4-10. Combined sine and cosine surface error. 


2 2 2 2 2 
Q-, = cos ip + sin <p - 2X^ sin (f cos cp 


( 4 - 47 ) 


Q ^ = (X + U )^ + (Y + V 
e ' s s s s 


( 4 - 48 ) 


The sagitta expressions are: 


= 






2 ^ 2 ^ 1/2 + + ... ( 4 - 49 ) 


S = 


I 2 
^2 ^2 




( 4 - 50 ) 


^2 - 1 - 


( 4 - 51 ) 


/ sin C ttQ \ 

^ = C, ( -pz Icos C_ ttQ 

' H *=2’'°a / ^ ' 


( 4 - 52 ) 


Note that Cj is an amplitude term while and control 
frequency. 

The ray- tracing equations for these special surfaces are: 


^s(i+l) 


^si- 


z . 

si 

1 - 




( 4 - 53 ) 
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^s(i+l) " 


where 


and 



2 

■ 

» 2 2.2 

1 +[;' (X)] 2 + [;' (Y)] 2 

1 9 S 


2 

"si 

1 + [C^(X)] ^ 

, (Y)] ^ 

-‘s0^s™-”s0^s 

(Y) 





(4-54) 


^'l<V 

(4-55) 

II 

£ 

*. w 


(4-56) 

(X) = p'l 

2 2 1 
X^^ (a + b sin (p) + b sin q) cos (p j 

(4-57) 

?.'l (Y) = Pi 

r 2 2 

1 (a + b cos (p) + b X^^ sin <p cos tp 

(4-58) 


n P'l 

(4-'59) 


. I 

"1 - Qj dQ, 

(4-60) 

^(X) 

” / 2 \ 
= Pt ( X . cos w - Y . sin (p cos (pj 

M \ S X S X / 

(4-61) 

^(Y) 

= P 2 (^si ^ " X^^ sin q> cos q^ 

(4-62) 


ft 


PZ = 





(4-63) 
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(4-64) 


I 



(Qa> = 


The refraction equations corresponding to eqs. 4-Z2 and 4-23 

are: 


I 


«si = ‘s-l 
""si 

or instead of eq 4- 22a and 4-23a 


(4-65) 

(4-66) 


£. = 2 - (X ) n . 

SI so S s SI 


(4-67) 


m . = m 
SI so 


C (Y )n . 
^ s s sx 


(4-68) 


d. Perfect Surfaces — Another kind of special surface 
which has proved of considerable utility is the option of interspersing 
one or more "perfect" surfaces in a system being analyzed. The most 
obvious application of this is the conversion of an afocal system into 
an image-forming system without introducing additional aberrations. 
The ray trace equation for this surface utilizes the standard intercept 
computations for a conventional plane surface. 


X 

s 


X 


s- 1 


f s- 1 

+ (t 


n 


s- 1 


1 - Z ,) 
s- 1 s- 1 


(4-69) 


Y 

s 


m , 

Y - + — ^ (t , - Z , ) 
s-1 n , s-l s-1 
s- 1 


(4-70) 


Z = 0 
s 
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The refraction equations are 


^S-1 

— = - P X 

"s *^s-l ® ® 


(4-71) 


m m , 

s _ s- I 

n n 


s- 1 


- P Y 
s s 


(4-72) 


n 

s 


jjs 

[l + (^g/n^)^ + (m^/Hg)^ 



(4-73) 


where P is the power of the perfect surface and it is input in place of 
s 

(RH(/JS) or as RADS=R^ = 1/P^. A perfect surface is identified by 
the input of the rotational symmetry constant (CRS) for that surface 
a s C RS = - 1 . 

e. Applications — Except for perfect surfaces, these equa- 
tions can represent most, if not all, optical surfaces generated by 
design or inadvertence. This includes spheres, cylinders, toroids, 
axicons, roofs, rotationally symmetrical surfaces with small cylindri- 
cal error, and numerous combinations of these as well as the periodic 
errors. 

The matrix in Table 4-1 illustrates some of the surface forms 
which may be represented by these equations. 

f. Optical Path Difference or Wavefront Phase Error — The 
optical path length D of a ray is the product of the length d along the 
ray and the refractive index p of the medium. 


D = optical path length 
D = d’p 


(4-74) 
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Defining the optical path length for a ray passing from surface s-1 to 

surface s to be D , it can be shown that 
- s 


D 

s 


(t , - Z , + Z ) 
s-1 s-1 s 


n 


s-1 


H 


2 

s-1 


(4-75) 


The total path length through the complete system is the s\im 
of all the surface-to-surface values and is defined to be 




D 

o 


n 



(4-76) 


which is in practice the summation from entrance pupil to exit pupil. 

is the distance along the ray from the object point to the entrance 
pupil for a finite object, i.e. ^ ^ 0. 


D 

o 


M- ^ 
’^o o 


(L '^ 4 0) 
n ' o 


(4-77) 


For infinite conjugates, is the distance from the incident plane wave 
to the entrance pupil. 


D 

o 


= m y 
O'' o 



0) 


(4-78) 


where 

y^ E entrance pupil ray intercept. 

Note that for the infinite conjugate axial case m^ = 0. , hence 
= 0. The difference between the total optical path length D* for 
an arbitrary ray and that for a reference ray is the optical path 
difference or wavefront departure W for the given ray. 
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(4-79) 


W = D* - D* . 

c 


W can be expressed as the sum of a set of W's at each 


surface 


N 

W = D+ VD-D 

O /. s CO 


= D - D + 

O CO ^ ^ 

1 


N 

I 


N 

I "c, 

1 


D ) 

cs 


N 

= Y (D - D ) 
Z. s cs 

0 


(4-80) 


or, defining 


W = D - D . 
s s cs 


(4-81) 


Then, for the ray corresponding to exit pupil coordinates (u, v) 
the overall optical path difference for that ray is W(u, v) and 


W(u 


N 

,v,= y 


w 


(4- 80a) 


s = 0 


Denoting the ray parameters for the arbitrary ray with the 
subscript s and those for the reference ray with the subscript cs and 
substituting (4-75) into (4-81), 


W = 


2 

^^s-1 


s n , n , 
s- 1 cs- 1 


n , (Z - Z i) -n 1 (Z - Z ,) 
cs-1 ' s s-1 s- 1 cs cs-1 


+ ‘s-l - ”s-l> 


(4-82) 
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Generalizing (4-82) to allow treatment of tilted and/or decentered 
surfaces 


W = W 


so 


2 

^"s-1 


n m , 
s- 1 cs- 1 


n' , (Z - Z -n' , (Z - Z 

cs-ls s s-l'cs cs 


(4-83) 


where 


W 


2 

^s-1 


so 


n , n , 
s- 1 cs- 1 


(n 


cs- 1 


■n ,)t ,-n ,Z ,Z 1 +n ,Z 

s-1 s-1 cs-1 cs-1 s-1 s- 1 cs-1 


(4-84) 


The primed and starred quantities of equation (4-83) are defined 
in Section 4,6 on tilt and decentration. represents the path dif- 

ference from surface "s-1" to the nominal tangent plane of surface "s". 
The remainder of eq. (4-83) represents the path difference from the 
nominal, untilted tangent plane of surface "s" to the tilted and/or 
decentered surface itself. 

g. Perfect Surface OPD — When using the perfect surface 

option, the optical path difference for the surface is computed as 
follows: 


W 

s 


W 

so 


p p 
s'^s 



+ Y 
s 

f^s + 



n 

s 



+ Y 


cs- 


+ n 


cs 


n 

cs 


(4-85) 
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4.4 


VIGNETTING AND OBSCURATION 


a. Vignetting — Vignetting must be accommodated and 

LiACOMA handles it in two ways; at the entrance pupil or on a surface- 
by- surface basis. 

If the nature of the vignetting at the entrance pupil is known, 
this vignetting may be described by the input of a value V^ = VZERO 
or V^ = VPI representing the fractional vignetting, respectively, in 
the upper or lower portion of the entrance pupil. A different set 
(Vq, V^ ) must be input for each field fraction K. The interpretation of 
the magnitude of the input value of V is as follows: 

V = 0 or 1. 0 no vignetting 
0 < V <1,0 normal vignetting 

The other means for handling vignetting can replace or augment 
the pupil vignetting. This is accomplished by the input of clear aperture 
data for each surface. As each ray is traced through the surface, a test 
determines whether the ray falls outside the clear dimensions in which 
case it is terminated. 

b. Obscuration — As with vignetting, there are two 

approaches to obscuration, either in the entrance pupil or on a surface - 
by- surface basis. At the entrance pupil, the descriptive data are: 
UZERO for the radius of the obscuration; GZERO representing the 
distance from the center of the entrance pupil to the upper edge of the 
obscuration; and GPI, the distance from the center of the entrance 
pupil to the lower edge of the obscuration. A set (UZERO, GZERO, 

GPI) is required for each field angle to be traced. A circular centered 
obscuration woxxld have 

GZERO = UZERO = - GPI 
Other cases can of course be accommodated. 
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In case (GZERO-GPI) > 2 • UZERO, such as in Fig. 4-11, no 
rays will be traced through the racetrack- shaped obscuration. If 
(GZERO-GPI) ^ 2 • (BETAO + UZERO), the effect is that of a rectangu- 
lar obscuration as per Fig. 4-12. 

When UZERO > (GZERO-GPI), the effect is that of a rectangular 
obscuration rotated 90° from that of Fig. 4-13. 

The other mode of obscuration input is via obscuration at each 
surface. 

c. Surface by Surface Vignetting and Obscuration — The 

parameters BETAPX, BETAPY, and BETASX and BETASY can be 
used to describe a circular obscuration or clear aperture, rectangular 
obscuration or clear aperture, annular obscuration, and segmented 
apertures. These are illustrated by Figures 4-14, 4-15 and 4-16. 

Some of the possible aperture configurations which can be described 
are shown in Fig. 4-17. 

4.5 REFERENCE WAVEFRONT 
th 

The N surface to which the OPD values are summed is 
analogous to the gaussian sphere discussed in the previous literature. 

The relevant characteristics of this reference wavefront or gaussian 
sphere are obtained from the chief ray data. As with any other optical 
surface, the reference wavefront is defined by its radius of curvature 
RW, its spacing TR from the (N-l)st surface, its decentration AX, 
ay and tilt 0X, 0Y. 

For the axial case where there is no tilt or decentration of the 
reference wavefront, its location is given by TR = Texit. 

Where Texit is due to eq. 4-6 and RW = T , - TR 

^ n-1 

Where T^ ^ was determined by eq. 4- 14a or eq. 4- 15a. 

Off-Axis Wavefront — For an optical system with tilts and/or 
decentrations or with an off-axis input beam the fvill set of descriptive 
parameters for the reference wavefront will be required. 
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Fig. 4-11. 

Racetrack obsuration (U = 
UZERO). 



Fig. 4-12. 

Rectangular obscuration 


3079 -3Z 



Fig. 4-13. 

Rectangular or near- 
rectangular obscuration 
oriented 90^ away from 
Fig. 4-12. (G = 
GZERO, G = °GPI). 
NAS 5-21898 ^ 

Final Technical Report 
Figs. 4-11, 4-12, and 
13. 
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OBSCURATION RADIUS = RQgg 
8ETASX = Rq^ and BETASY = 0 

OR BETASY = R^,^ AND BETASX =0 
BETAPX = Rggg AND BETAPY= 0 
OR BETAPY = Rq35 AND BETAPX= 0 


Fig. 4-14. Circular clear aperture and obscuration. 
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Final Technical Report 

Fig. 4-14 
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BETAPX = Xqbs 
BETAPY = Yqbs 


Fig. 4-15. Rectangular clear aperture and obscuration. 
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Fig. 4-15 
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ANNULAR OBSCURATION 



RAYS TRACED INSIDE R|^,|^ OR OUTSIDE R|viax 
RAYS BLOCKED BETWEEN Rn-,|>, AND R^ , y 
BETAPX = R|y,„y, AND BETAPY = 0 

OR BETAPY = R|y„fg AND BETAPX = 0 
BETASX = AND BETASY = 0 

OR BETASY = Rj^^^ BETASX = 0 
(ONE OF THE ABOVE VALUES MUST BE INPUT NEGATIVE) 

Fig. 4-16. Annular obscuration. 
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Fig. 4-16 
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02 



RW is the distance along the chief ray from its intersection or 
point of closest approach to the optical axis (exit pupil position) to its 
point of intersection with the image surface. The angular components 
of the chief ray give the tilt of the reference wavefront, etc. The com- 
putation of RW, T^, 0X, 0Y, AX and AY is as follows. 

Given 2^, m^, n^, for the chief ray exiting the optical train and 

the chief ray intercepts X , Y and Z at the focal surface, if m 4 0 
^ ^ c c c c 


RW 



(4-86) 


TR 


T • 
n- 


1 


Y n 


c c 
m 


c 


+ Z 

c 


(4-87) 


^ AX 


X 

c 


fcjc 

m 

c 


AY = 0. 

If m =0 and 2 ^ 0 

c c 


(4-88) 


RW 


X 


n- 1 



(4-89) 


X n 

TR = T' , ^ + Z 

n-1 c 


(4-90) 


AX = 0. 

m RW 

AY = Y ^ 

^"s-l 


(4-91) 
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The reference wavefront tilt angles are 


ex = ARCTAN (i /n ) 
c c 


(4-92) 


0Y = ARCTAN 


m 


L"c(‘ 


(4-93) 


For a case with an IMC so that the reference wavefront is off- 

axis but X — Y — * 0, TR is defined as that for the axial case and 
c c 


RW = 


(T' , - TR) 

n- 1 


Hn-1 


(4-94) 


AX = X 


H RW 


(4-95) 


AY = Y 


m RW 
c 

^N- 1 


(4-96) 


Equations 0X and 0Y are defined per eqs. 4-92 and 4-93, 

4.6 TILT AND DECENTRATION 

The tilt and decentration procedure utilizes conventional 
coordinate transformation matrices. 

The procedure is to trace a ray to the tangent plane for the cen- 
tered surface and then locate the intersection point of the ray with the 
tilted- dec entered tangent plane, transform the optical direction cosines 
of the ray into the rotated coordinate system, and then trace the ray 
from the new tangent plane to the surface. 
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The transformation matrix represents the relocation of the 
origin by AX, AY and AZ followed by rotation of this displaced coordi- 
nate system. The rotation begins with a 6Z around the optical axis fol- 
lowed by the rotation 0Y around the X' axis where the X' indicates the 
new X axis due to the Z rotation. 0X is the rotation about the Y' axis 
which is the result of the 0Z, 0Y rotations. 

Representing the new optical direction cosines with primes, 

O' I = 2 1 (cos 0X cos 0Z + sin 0X sin 0Y sin 0Z) 

'^s-l s-l' 

-m , cos 0Y sin 0Z + n , 
s- 1 s- 1 

(cos 0X sin 0Y sin 0Z - sin 0X cos 0Z) (4-97) 

m ' , = m . cos 0Y cos 0Z - n , (sin 0X sin 0Z + cos 0X sin 0Y cos 0Z) 

s-l s-1 s-1' 

+ 2 , (cos 0X sin 0Z - sin 0X sin 0Y cos 0Z) (4-98) 

s- 1 

n' , = n , cos 0X cos QY + I , sin 0X cos 0Y + m , sin 0Y . 
s-1 s-1 s-1 s-1 

(4-99) 

The values X=I% Y'*' and represent the intercept coordinates 
of the ray with the nominal tangent plane expressed in the transformed 
coordinate system relative to the new tangent plane. 

X^' = (XSO-AX) (cos 0X cos 0Z + sin 0X sin 0Y sin 0Z) 

-(YSO-AY) • cos 0Y sin 0Z - AZ 

(cos 0X sin 0Y sin 0Z - sin 0X cos 0Z) (4-100) 
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Y* = (YSO-AY) cos 6Y cos 6Z +AZ (sin 6X sin 6Z + cos 6X sin 0Y cos 0Z) 
+ (XSO-AX) (cos 0X sin 0Z = sin 0X sin 0Y cos 0Z) (4-101) 

Z>!= = (XSO-AX) sin 0X cos 0Y + CYSO-AY) sin 0Y - AZ cos 0X cos 0Y . 

(4-102) 

Thus X' and Y' represent the intersection of the ray with the new 
so so ’ 

tangent plane and 

, ^s-1 

X = X» - f ■ (4-103) 

so n* , ' ' 

s- 1 


Yt 

so 



(4-104) 



a. Sign Convention — The origin (surface vertex) of the trans- 

formed coordinate system has the coordinates AX, AY, AZ relative to 
the original origin (Fig. 4-18). Viewed along the positive direction of 
the optical axis, +0Z represents a counterclockwise rotation about the 
optical (Z) axis. Similarly, a tilted surface has its optical axis oriented 
at the angles 0X and 0Y with respect to the axis of the preceding optical 
surface. A surface tilted by +0Y (or + 0X) has it axis rotated such that 
its new direction relative to the original axis is the same as would a 
ray having positive direction cosines m and n. For example, a ray 
parallel to the original axis (£ = m = 0, n = 1), and incident upon a 
surface tilted by a positive 0Y will have transformed incident direction 
cosines so that m < 0 (see Fig. 4-19). 
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Fig. 4-18. Illustration of AY negative. 
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Fig. 4-19. Illustration of 0 positive. 


NAS 5-21898 
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Figs. 4-18 and 4-19 
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b. 


Coordinate Restoration — If a surface is tilted or 


decentered, all the surfaces following will be aligned with the tilted or 
decentered surface. Therefore, the tilt or decentration of a single 
surface amounts to tilting or decentering the system following the tilt 
with respect to the portion of the system preceding the tilt. If it is 
desired to misalign a given surface in the midst of a centered optical 
system, it is necessary to introduce a reverse tilt or decentration 
immediately following the misaligned surface to restore the coordinate 
system. If there is a combination of AX, AY, AZ and/or 6X, 0Y, 0Z 
at the surface, then because of the sequence of the operation of the 
decentrations and tilts the combination of - AX, -AY, - 0X, -0Y, etc. 
will not necessarily restore the proper coordinates. The following 
are rules to follow in the restoration of the coordinate system after 
traversing a misaligned surface: 

1. Restoration after a single tilt 0Y or a single 
decentration AX is accomplished with the inser- 
tion of a dummy surface following the misaligned 
surface with = 0, Tg_ j = 0 and the refractive 
index Pg the same as that following the misaligned 
surface. The dummy plane should have a tilt of 
-0Y or a decentration of -AX, 

2. For a surface having tilts 0X, 0Y, 0Z and decen- 
trations AX, AY, AZ. Four dummy planes must 
be inserted: 

Dummy plane No. 1 with -0X 
Dummy plane No. 2 with -0Y 
Dummy plane No. 3 with -0Z 
Dummy plane No. 4 with -AX, -AY, -AZ 

3. For surfaces with more than one but less than the 
full six misalignment components the sequence of 
the operations must remain the same removing 
dummy planes which do not apply. For a surface 
with 0X, 0Z and AY three dummy planes are used: 

Dummy plane No. 1 with - 0X 
Dummy plane No. 2 with - 0Z 
Dummy plane No. 3 with -AY 
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PUPIL, COMPUTATIONS 


4. 7 


All the ray trace and OPD computations are initiated at the 
entrance pupil. Part of the input is "M" which is the surface n\amber 
of the aperture stop which provides for the computation of the paraxial 
entrance pupil location T^ from eq. (3-5). The rays are input at the 
pupil in a uniform square grid as in Fig. 4-20. The input parameter 
NXY determines the number of grid intervals across half the pupil so 
that there will be (2NXY + 1) grid points across the pupil. The pupil 
semidiameter is divided by NXY to give 


ex = ey = 


pQ 

NXY 


(4-105) 


for the intervals between rays. The ray coordinates X and Y at the 

entrance pupil are then obtained by incremental steps of ex and ey. 

Combining these with the input values of L^ ^ and H for any angular 

obliquity of the input beam the ray input direction cosines 1^, m^ and 

n are determined 
o 


AA = ARC TAN (X L '^) 

o o 


(4-106) 


BB = ARCTAN ((Y^ L^'^ + Tan (H)) cos (AA)) (4-107) 


1^ = sin (AA) cos (BB) 


(4-108) 


m^ = sin (BB) 


(4-109) 


n^ = cos (AA) cos (BB) . 


(4-110) 


At this point is determined to be used in equation (4- 80a) 
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Fig. 4-20. Uniform square ray grid. 
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(4-111) 


(4-112) 


Pupil Amplitude Values (B(u, v)) — The pupil coordinates are 
here referred to as u, v which coincides with the notation used in dis- 
cussing computation of the amplitude spread function (Section 4, 8 ) u 
and V apply to the exit pupil coordinates while and represent the 
entrance pupil coordinates for the given ray. 

The entrance pupil amplitude values are determined as follows: 

Transmitter — Truncated gaussian amplitude distribution with 
total radiant power per input parameter P^. 

Receiver — (Received signal) — Uniform amplitude distribution 
if Pq is input positive it is the uniform intensity value If P^ is input 
negative it is the total radiant power, then 

I = Ip 1 /A (4-113) 

o ' o' o 


where 

= entrance pupil area 


then 

B(u, v) = . 


Receiver- Local Oscillator — Truncated gaussian amplitude dis- 
tribution with total radiant power per input parameter P^. 


Ill 



Truncated Gaussian Intensity Function — The pupil amplitude 
function B(u, v) for a gaussian amplitude distribution is defined as 
follows: 

= Total input beam power (W) 

= peak intensity of gaussian beam 
2 

= l/e intensity radius 


= maximvim input beam radius or truncation radius 


2 2 
r = u 


+ V 


2 


I 

o 


2 

O 


(1 


2P 

o 




(4-115) 


I(u,v) 


= I e 


o 



(4-116) 


B(u, V,) = yjl (u, v) . 


(4-117) 


The amplitude expressions, eq. (4-114) and (4-117), assume unit 
magnification between entrance and exit pupils. If the pupil magnifica- 
tion mP is not unity it will modify B(u, v). 


pupil magnification = mP = 
b'(u, v) = mP 


entrance pupil diameter 
exit pupil diameter 

• B(u, v) . 


(4-118) 

(4-119) 


4. 8 AMPLITUDE SPREAD FUNCTION 


The complex impulse response function or amplitude spread 
function ASF(x, y) at the Fraunhofer plane is the two dimensional Fourier 
transform of the pupil fxinction F(u, v). 
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For the received signal 


ASFl(x,y) = FRl(x,y) + i FIl(x,y) (4-123) 

PSFl(x,y) = |ASFl(x,y)|^ (4-124) 

and for the local oscillator 

ASF2(x,y) = FR2(x,y) + i FI2(x,y) (4- 123a) 

PSF2(x,y) = |ASF2(x,y)l^ . (4- 124a) 


Signal power, local oscillator power and heterodyne power are 
obtained by numerical integration over the specified detector dimensions. 

The computed ASF and PSF values are obtained for a uniform 
grid of intervals AX, AY. 

Signal Power = POWER(l) = Z S PSFl (x, y) AX AY (4-125) 

X y 

L.O. Power = POWER (2) = ZZPSF2 (x,y)AXAY . (4-126) 

X y 

Heterodyne Power — If the nominal phase difference between the 
receiver signal beam and the local oscillator beam is zero, heterodyne 
power is obtained with J12 (x, y). 

Heterodyne Power = POWER(3) = Z Z J12(x, y) AX AY (4-127) 

X y 


where 


J12(x,y) = FRl(x,y) FR2(x, y) + FIl(x,y) F12(x,y) (4-128) 
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ASF(x, y) = F(u, v) 


(4-120) 


^denotes Fourier transform 

F(u, v) = B(u,v) • exp [-iK • W(u,v)] (4-121) 

K = 2n/k . 

With B(u,v) from (4-114) or (4-117) and (4-119) with W(u,v) 
from (4- 80a). 

ASF(x,y) = — ■ exp |-^|^(ux + vy) j dudv. (4- 120a) 

The complex amplitude ftinction ASF(x, y) is represented by the 
real function FR(x, y) and the imaginary function FI(x,y). 

ASF(x, y) = FR(x,y) + i FI(x, y) . (4- 120b) 

The intensity function or point spread function is PSF(x, y) 

PSF(x.y) = |ASF(x,y)|^ (4-122) 

= [FR(x,y) + i FI(x,y)] [FR(x,y) -iFI(x,y)] . (4- 122a) 

The functions ASF(x, y) and PSF(x, y) provide the basis for the 
quality criteria for the receiver and transmitter. 

4.9 RECEIVER QUALITY CRITERIA 

The cimplitude spread function ASF(x, y) and point (intensity) 
spread function PSF(x, y) are determined for the received signal beam 
and the local oscillator (L.O. ) beam. 
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or for brevity, 


J12 = FRl FR2 + FIl F12 . (4- 128a) 

If there is a nominal phase difference (A<j)) between the two beams, that 
phase component must be removed to obtain J12: 


Heterodyne power 


POWER(3) =22 J12' (x, y) AX AY 
X y 


(4-129) 


J12' = FRl' FR2 + Fir F12 (4-130) 

J12' = (FRl cos A<j> - FIl sin A<t>) FR2 
+ (FIl cos A<^> + FRl sin A4>) FI2 

J12' = (FRl FR2 + FIl FI2) cos A<j> 

+ (FI2 FRl - FIl FR2) sin A(j> 

J12' = J12 cos A<i> + K12 sin A<)> (4- 130a) 

where 

K12 = FI2 FRl - FIl FR2 (4-131) 


Phase minimization — The value of J12' will be maximized when 
d(J12'/d(A(j)) = 0. 


d(J12') _ 
d (A4>) 


J12 sin A<|> + K12 cos A<t> 


(4-132) 


with 


d(J12') 
d (M) 


0 . 
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J12 sin A4> = K12 cos A4> 


sin A(t > _ K12 
cos A(j> ~ J 12 ~ 


A4> . 


The value of A<f) is then obtained by the summation over the 
specified detector as follows 


Tan = 2 2K12 (X Y) 

X v 

^ ^ J12 (X Y) 
X y ' 


(4-133) 


The remaining received signal quality criteria are determined 
as follows: 


Phase Match Efficiency 


POWER (3) 
POWER (5) 


(4-134) 


POWER (3) = Heterodyne power 

POWER (5) =ZZ [PSF1(X, Y) PSF2(X, Y)j 
X y 


1/2 


AX AY 
(4-135) 


Optical transmission 


POWER(6) 

PZERO 


P0WER(6) = Exit pupil power 
PZERO = Total received power 


Focusing Efficiency 


POWER(l) 

POWER(6) 


(4-136) 


(4-137) 


Nonheterodyne detection efficiency 


POWER(l) 

PZERO 


(4-138) 


L.O. illTimination efficiency 


POWER(5) 

(POWER(l) POWER(2)^^^ 


(4-139) 
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(4- 140) 


McLxirmim Antenna Gain 


4tt 

o 



A^ = entrance pupil area 
X. = spectral wavelength 


Receiver efficiency to I. F. 


(POWER (3))^ 
PZERO -POWER(2) ’ 


(4- 141) 


4. 10 TRANSMITTER QUAUTY CRITERIA 


The point (intensity) spread function PSF(X, Y) is determined for 
the transmitter at the specified range. The quality criteria for the 
transmitter are: 

Peak intensity = IP = maximum value of PSF (X, Y) (4-142) 


Optical transmission 


POWER( 6 ) 

PZERO 


(4- 143) 


Maximum antenna gain = — POWER(7) ^ 44 ) 


POWER(7) = Effective exit pupil area 
= Spectral wavelength 


Overall transmitter efficiency 


IP (RW xf- 
PZERO POWER(7)' 


(4- 145) 


where RW = Range. 
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4.11 ERROR ANALYSES 


The errors associated with these computations come from 
numerous sources. The major error sources are discussed here with 
some indication of guidelines to be followed to keep them to a minimum. 

a. Aberrations — The conventional Fraunhofer diffraction 
integrals used in these computations are derived based on certain 
approximations. The f\indamental integral involves the exponential 

exp [iK^W + S)] (4-146) 

where W = OPD error for a given point on the wavefront and S is 
defined to be the distance from the gaussian sphere to a given observa- 
tion point in the intensity (or amplitude) plane as illustrated in 
Fig. 4-21. 

The assTxmption is that S *R in the denominator of eq. (4-146) 
where R = radius of gaussian sphere. It is also assximed that in the 
numerator of eq. (4-146) 

S.R - 

where u, v represent the coordinates of a given wavefront point and x, y 
are the coordinates of the observation point in the intensity plane. As 
the aberrations become large, these approximate representations begin 
to break down and the computed results lose accuracy. There is no 
well defined threshold of aberrations where the diffraction integrals 
become inaccurate; but, in general, it can be said that results are best 
when W/X. « 1. 

b. Aperture Description — The pupil amplitude B(u, v) and 
OPD W(u, v) data are obtained for a rectangular array of points uniformly 
distributed over the pupil. For a circular aperture the process of fitting 
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Fig. 4-21. Deviations from the gaussian sphere. 
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square pegs in the rovind aperture hole can only represent the aperture 
botindary accurately if there are a large number of small square pegs. 
This effect is further compoxinded by circular obscuring apertures, or 
even rectilinear blockages such as spiders or support struts which do 
not line up with the u, v coordinate system. Figure 4-22 illustrates how 
the rectangular u, v array describes a circular aperture boundary. The 
uncertainty of determination of the area of the pupil, hence signal 
power, etc. will be ± dA where 

dA = tX eY (4-147) 

There are techniques for weighting the amplitude values of array ele- 
ments lying on the boundary for which there is an accuracy improvement. 
The tradeoff of accuracy gain versus increased computational complex- 
ity of these weighting techniques must be compared against a simple 
increase in grid size (number of array points). 

In terms of some representative values of NXY, the normalized 
uncertainty dA/A woxild vary per following tabulation. 


NXY 

dA/A 

10 

0. 0032 

15 

0. 0014 

20 

0.0008 

25 

0. 0005 

30 

0. 0004 

40 

0. 0002 

50 

0. 0001 


c. Circxilar Detector — The errors associated with a circu- 

lar detector sampled via the square spread function grid will follow the 
same rules as the aperture description. The \incer tain ties in the com- 
puted res\ilts associated with the detector would follow Table 4-2 
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with NXY replaced by IDIM. IDIM is the array size for the 
detector. 


IDIM = 


Circular Detector Diameter 
2 DELFS 


IDIM = IDETEC (2, J) 
= 10 (Default) 


d. Aliasing — The computed Fraunhofer diffraction function 
is based on a discrete complex amplitude array (pupil function). This 
discreteness results in a series of replicated, overlapping far-field 
patterns as illustrated in Fig. 4-23 which constitute an infinite periodic 
f\anction. The spacing between these replicated patterns is directly 
proportional to the array size or number of points across the pupil. 
These patterns overlap so that each computed ASF or PSF value repre- 
sents the superposition of all the replicated spread functions at a given, 
point in the image plane. Thus, a point midway between two replica- 
ted spread fvinction peaks will have a value twice that of either spread 
function plus a contribution from all the other surrounding spread 
functions. The larger the number of sample points across the beam, 
the greater the separation between these replicated spread functions, 
hence, the lower the errors in the spread function due to aliasing or 
overlapping. 

Aliasing errors increase with aberration since the effect of 
aberrations is to increase the spread function values away from the 
center while reducing the central spread function values. 

e. Effect of Defocusing on Aliasing — As indicated in 
Fig. 4-23, the distance between peaks of the replicated overlapping 
spread fiinctions is Axrep and 

Axrep = 2X NXY f# 


I 


I 
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Fig. 4-23. Replicated (aliased) spread functions due to discrete sampl- 
ing of pupil. 
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where 


X = spectral wavelength 

NXY E number of grid points over entrance pupil semidiameter 
f# = relative aperture 

Aberrations and/or defocusing which cause geometrical spread 
approaching Axrep will destroy the accuracy of these computations. 

A defocusing AF res\xlting in a geometrical spread of Axrep wovild be 

AF = 2 f# Axrep 
= 4X NXY (f#)^ 

which is obviously excessive as would be a defocusing of half or one- 
fourth that amount. It can be said, therefore, that to minimize aliasing 
error due to defocusing 


AF« X NXY(f#)^ 


and preferably 


AF — 0. 


f. LACOMA Test Case — Numerical Results — A test case 

was devised to check accuracy. It consists of diffraction limited 
received signal and L. O. optics. The relative aperture of the received 
signal optics is 10, and for the Li. O. is 100. The matrix of optical 
path differences indicates zero (10 over the pupil. The (complex 
amplitude) pupil function indicates a uniform 100%. A peak value of 
1. 0013013 rather than 1. 0 indicates a slight demagnification of the 
exit pupil as used for the analysis. The ASF for the received signal is 
real with a peak value of 73, 74928 compared to a predicted value of 


74. 0147 


a 

o 


A P 
c o 

^2 2 


1/2 


74. 0147 
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= entrajice pupil area 
= total input power 
F = focal length 

representing an error of 0.3%. Signal power is computed to be 
65.6453 which compares with a predicted diffraction limited value of 
65.82 for an error of 0. 3%. Both these errors are to be expected from 
the value of IDIM = 10 and the predicted uncertainty of 0. 0032 for 
NXY = 10 per Table 4-2. 

The local oscillator results were examined in terms of the peak 

ASF value a^. The computed value is 0. 8018 while the predicted value 

is a^ = 0.8352. The 4% error here is due to the gaussian pupil function 

changing rapidly at the edge of the aperture where the aperture sampling 

2 

errors are thereby enhanced. The further the l/e point is from the 
aperture edge or truncation point (in either direction) the lower this 
error will be. 

Review of numerical results for the various sample cases indi- 
cates that the above accuracy values represent a realistic range to be 
encountered in the analysis of typical systems. 
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APPENDIX A 


TRANSMlTTi^ TfST CASE ••• 

$INP 

B£TA0(1) s3., .3, K(ll s 3«., H(l) su.« 0* » VLlNVIlt s0.*3 
HAv(l) sC.iO.t IFLGKII s wtU« 

RH0S(1) 

TS(1> sc., 

XJIUS(l) = 0*1. , 

Ci?S(2) =-l., CJ?S(5) s -1,, 

NFRS = li], OELFS » .U, 

MXY 3 25, 

SEND 


I Preceding page blank 
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LACCMA 


LftSFR CDHMUNICATOR ANALYSIS PROGRAH 


PAGOS CORPORATION 


••• TRANSMlTTtR TEST CASE ••• 


N 


Bc-TAO S.UOwJC..JE«’UO H Q. 


P7ERO 0. 


M : . XLlNV 


HAG 


IFLGl 0 


SURFACE 

PARAMETERS 





SURF. 

CURVATURE SEPARATION 

INDEX 

OBSCURATION 

DATA 

NO. 

R«OS TS 

XHUS 

3ETAPX 


BETAPY 

i 

G. 0. 

I.COQjOJ 

0. 

C. 


2 

l.LwuwtiJt.E~w2 l.UUWJytlUE^UZ 

1. C30003 

J. 

C. 


3 

C. 0. 

l.CCObOO 

3. 

0. 



SPECIAL SURFACES 


VIGNETTING OATA 
BETASX BETASY 

b. 0. 

u. . 0. 

0 * 0 • 


SURFACE RH3i 


RH02 


CRS CC PHI 


AHP 


2 


-i.cco jooo o«oooc:cc 


0. OOOODC. 


C.JGOOOOO 0 


{ANGLES INFUT GREATER THAN 6.2£3 ARE ASSUHEO TO 8£ IN DEGREES) 


FIELD ANGLE « 0« 


NXY s 2S> NOET s 1 NFRS s 1) 

OELFS s !• uC03oO<)Lc*C2 ZOHEGA s C« VLAHOA s : 6ll3aSwe*02 

OETECTOR NO. 1 lOETEC * 1 0 0 0 0 


FREQ 

.0300000 
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P^RaXIAL PAY TRACE 


BETA 

R 

ALPHA 

A 

5*. COkClgOE^ul 

0. 

0. 

2.C JOOOOOCE*31 

5«.Lwclu03E*'-i( 

•5. o:Ogg jC 0£*02 

a. 

2.0G0O0;0UE-01 

G • 

*S«u03g0iLgE*02 

2.000C J J-jJE^Ol 

2.0 J00Q300E-31 


PARAXIAL entrance PfPiL POSITION T(0» = -0. TEXIT® 0. 

INVERSE OBJECT DISTANCE XLlNV = 0. 

FOCAL LENGTH s 1 . L J Oc J Jg 02 RACK FOCAL LENGTH s 1. OOC : OC) CE» 02 

TS(N-l) = 1 • Cg 3 j UjC iE^Q2 Fl = I* 00 0 OuQ J(i£^02 
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* TPflfiS'll TTc‘'< T£?T CAS£ 


• « * 


CHirp PA> UtTA F0*< TRANS'II ritR. 


FIELD ANGL£ <JAOIA^S> = J.oidgOuJL 
CHIPF RAY CLCPJIMATF5 ON SJ^FAC- N 

> Y 7 

0 . : . C • 

DIRECTION CC^JNES or. 5 ?urfa:£ N-1 

LIN 

0 . .. i*cjo£lju:e «^07 

RADIUS RH AnO position TR OF £*IT O'jciL 
i-M • TR 

l*C-L...buE^U^ j » 

tilt and displacement of iKlT P'JPIL IN y- ANO Y-OIRECTIONS 


XTILT 

c. 


YTIlT 


XOISP 


YOISP 

3. 


KAVfcFROKT ERROR INFOK-IATION 


VARIANCE » U.iw.23327c-l-» 

RMS = **. 3C3?u957E-1h 

NAV.ELENGTH = l.^6U385St-J2 

MAKIMUH UNNORMALl2tr ERROR* i . 33226763E- 1 5 
MINIMUM UNNURMALIZEC ERROR «- I . 33226763E* 1 5 
MAXIMUM NORMALI2EO ERROR* 1 . 255 5 . 7 776- I 3 
MINIMUM normalized 6 RPOR*- I ♦ 255 5^77 7t- 1 3 
APPROXIMATE STREHL PATIO * i . d 0 d 00 3 J 3£ ♦ C 0 
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TRft-JSMITTciR TiST CAS- ••• 

TRANSHITTER OATA. 

PEAK INTENSITT = t . 41223 19£ ♦ w I 

optical transmssion s 9 . j-iSbeabe-Qi ( ••C02 oe i 

NAXIHUH AnT'NNA GAtr 3 d.^3ll969£*i6 ( 69.411 OB > 

OVERALL TRANSMTTSR cPPlCIiNCY s 9. 22df 7 33E«tl t -.349 OB ) 
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APPENDIX B 


t«iansmitter test case *l 

SINP 8ETAl(1) s 2.23* HUI s j., IFLCKU s Q, n<1I s -7, 
RAOS(t> « I.., (i3«l<*952, o7.b95^4, 0 . • • 1»; 2d. , 6 .0£^6 « J.« 
TS(l) 3 2.«3e.ll« -Sm.t d«, 0.* 

20MEGA = d.925, 

XrfUS(l) 9 1., ^.0di62t l.f i.f -l.t -t.» 

SKAPA(5» 5 J,, 

CRS(6» s -i. , 

MFRS = 51| OELFS = 50t 
lENO 


I PreGeding page blank 

I 
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• • ♦ ♦ LCCOHA 


lAScR CJrlMUUlCarOR analysis procrah 


PA«05 C0RP3RATI0N 


T lANSMlTTfc:;< TISl CA$c *i 


N 


-7 


R* TAfl 


2«?90Cu(jJj.*J0 


H 0. 


P7CRO 0. 


H 0 XLlNV Xi* 


MAG 


IFLGl 0 


SURFACt PARAMETtRS 


SURF. 

CURVATURE 

SEPARATION 

INDEX 

OBSCURATION DATA 

vignetting data 

NO. 

RHCS 

rs 

XMUS 

RETAPX 

BETAPY 

6ETASX 

BETASY 

1 

c. 

4*14jjDw3E^02 

1.CQ0..03 

J. 

0. 

c. 

0. 

2 

2s3l/5229E- J2 

2«auai4uuE^^(i 

4.G3362) 

3. 

c. 

0. 

0. 

3 

1.4/du77>£-Ci 

3.611job JE^Ol 

l.COObOQ 

0. 

c. 

c. 

0. 

4 

C. 

S.14jJjo3E4b2 

I.OOOlQU 

U . 

c. 

c. 

0. 

5 - 

9. 7276265E-24- 

5. 14jJ3uoE^b2 

-i.coaooo 

tf . 

c. 

- G* 

0. 

6 

1.25.wuujE-'j7 

0. 

-i.coouo: 

u. 

c. 

0. 

0. 

7 

0. 

0. 

•l.C DOJOT 

0. 

c. 

0* 

0* 

SPECIAL 

SURFACES 







SURFACE 

PHOl 

RH02 


CRS 

cc 

PMI 

AMP 


FREQ 


6 


Ci. j. -i.&oooouo (l•oo!)o;ob o.oiooo:. 

(ANGLES INPUT GREATER THAN 6.263 ARE ASSUMED TO PE IN DEGREES) 


C.OCOOdOO 0*0000000 


SURFACE SFAPA 


-——ASPMERIC COEFFICI ENTS- 

ALPHA SETA GAHHA DELTA EPSILON 


6 


0. 


field angle 8 3. 


NXY * 2S NDET * 1 NFRS = 51 

OELFS s 5. :COyUuS(£vQi ZOMEGA « 9* 25031 OOOC-01 VLAHDA* 1 » ^ 61i365bE-0 2 
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PaRAXtAL PAf TRACE 


PETA 

2 . 25 wt.jOo£ 4 -Jl 

Ca2^W^t}GJt}£^uC 

2, 17177956£*CI 
-‘..4i*2e2776£-07 
-3.v9lj734lE»jl 
-3.i;9l27337£*-01 
-3.09l57337£fjl 


P 

3 . 

-i.56465l24e-0i 
-6.0l4345C7£-02 
-6.0143451 re-02 
d.64362aLHE-lQ 
-3.66335235E-06 
- 3 .a 6335237 e-Q 6 


ALPHA 

C . 4 

l.d4QCj0C0£«‘C2 -1 

1.77a25494£«>02 -4 

1.6626661 3£*>01 -4 

-2.27791996£*^C3 -3 

-2.294546a6£*>03 -3 

-2.2945-«&a6£«^C3 -3 


PARAXIAL ENTRANCe PL'PIL POSITION T(0l 


TEXIT= 7,: 3u97397£^34 


iNVtRSt 06JLCT QISTANCt «LtMW * 0. 


FOCAL LENGTH s 5 . 6239575 3E^ 35 BACK FOCAL LENGTH s 6. 00 1 7 S967efQ6 

TS(N-l) = b.LJl78967e>06 FL = 5 . 62 395753£^35 


A 

t*4444444£-3 1 
23509257£^01 
464C9694e»00 
46409694E>30 
234aCiiaiE-02 
26346365E-32 
26346365E-02 
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T-lANShlTTE? TEST CASE *1 


CMIEP RAY DATA FOR TRANSilTTER. 


FIELO ANGLcI (RAOIAKS) = u.JiJljOJlO 
CHIEF RAY COORDINATES ON SJRFACZ N 

X Y 2 

0 « * • C * 

DIRECTION COSINES ON SURFACE N- 1 

L i N 

C* J. •l.O'JCOODOOE^C) 

RADIUS RH AnU POSITION TR OF £xIT PUPIL 

TR 

7,9U'-*6C13£*36 7 . Q 3 397397£4 O*. 

tilt ANO OISPLACtMcNT OF £<XT PUPIL IN Y- ANO Y-OIRECTIONS 


XTILT 

C. 


ttilt 


XOISP 

C. 


YOISP 

3. 



wavefront error Information 


variance = H.97S3627LE-J4 

RMS s 7.Ad377j05E*GL 
WAVELENGTH = 1 . 361 13 8S0E- ^ 2 

MAXIHUH UNNCRHALlZEf- ERROR* 1 . 3 61122 76£*1 0 
MINIMUM UNNORMALI2EC ERROR*-l . 77fl4*6695c-fl5 
MAXIMUM NORMALIZED ERROR* 1. 301SL6 jaE* 16 
MINIMUM NORHALIZEO ERROR**!. 666<«6l7lt*a3 
APPROXIMATE STK£HL RATIO * 9. 99977889E-C1 
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TRANSMITTER TEST CASE *l 


TRANSMITTER OATA. 

PEAK INTENSITY = 1 . Aid? JJUE-Q T 

OPTICAL TRANSMISSION 3 9. ISSSSSAE-0 1 ( -.CllO 09 t 

MAXIMUM ANTENNA GAIK = 3. 3319166E*a8 ( 9S.2E? 09 I 

overall transmitter efficiency a 3.366C??4E>Ct « -4.729 09 ) 
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APPENDIX C 



RECEIVER TEST CASE *l 

tINP PcTAi,<ll 3 d2«9tl*39, H<t» s Q.«C.» XLiNVfU « u*fO.» »<AG(i) s 
N(l) a -9,-7* tUl\ a 5,6« IPLCKll a l,a* 

RAOSdl a w.t -320.«0.« l4l.aQi75t 6«0.» 0.« 

TSdl a 3b6.723« G«« •5:2.9J983, 1(*6. 1 A 3 S7 , 96 .423 • 2.t -2« • 

TS(6) - 6b«773« u*t Q** lll*2t G«« 

XMUSd) a l.t l.» i.i 4.JU062» 6*t«* 

SKAPA<2) a a., 1., .7173b, 

CRS<11) a -I*, 

TMTAY(b) a ’*.78539816, b«, -.78539816, 

BETAPxd) a 12.45, 

M<1) a 5,1 l, 

IPROGdl a 1,1, l,J, 

SEND 


Preceding page blank 
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• ♦ • • • LACOHA 


LASER COMMUNICATOR ANALYSIS PROGRAM 


PAGOS CORPORATION 


RECEIVER TEST CASE *1 


RECEIVEO signal OPTICAL OATA. 


N 

-9 

BETAO tt.290u'J.^u£«-0i 

M 3, 

P2ERO 0. 

H 

5 

XLINV C. 

MAG 0. 

IFLGl 1 


SURFACE parameters 







SURF. 

CURVATURE SEPARATION 

INDEX 

OBSCURATION DATA 


VIGNETTIN 

G DATA 

NO. 

RHOS T5 

XMUS 

RETAOX 

SETAPV 


8ETASX 

SETASY 

1 

C. 3.b6725JCE«^02 

I.COOjOO 

1.2453E*C1 

0. 

C. 


0. 

2 

-1.2195122E-J3 0. 

•1.G30J00 

w • 

0. 

G. 


0. 

3 

0. -5.329J9A9E«32 

-1. LOO. 33 

a. 

c« 

C. 


0. 

4 

7.G52..99iE*Q3 1 . ••614 3 A7E *^C2 

1. CGoaoa 

3. 

c. 

0. 


0. 

5 

C. 9.a**23JU0E«>01 

1.C00CI3Q 

j. 

c. 

s • 


0. 

6 

0. 2 • <iG«i J J u GE^CO 

4.C00&23 

3. 

c. 

c. 


0. 

7 

0. -2.0C3g00O£«-GC 

I.COOjGO 

3. 

c. 

0. 


0. 

6 

G. S.6773a00E«ul 

1.C3QOOO 

J. 

c. 

*■ 


0. 

9 

C. 0. 

L.lG3»0d 

0. 

G. 

w • 


0. 

TILT 

AND CISPLACEHENT DATA 







SURF. 

THTAX TMTAY 


THTA2 

DLTAX 


DLTAY 

DLTAZ 

6 

L. 7.8539616E* 

0. 

C. 


c* 


0. 

6 

u, •7.d5398l6C« 

ill a. 

0* 


3. 


0. 


(ANGLES INPUT greater THAN 

6.263 ARE 

ASSUMED TO BE 

IN OtiREESl 












surface skapa alpha 


BETA 

GAMMA 


DELTA 

E*»SILON 

2 

C. u. 

0. 

0. 


.0. 


0. 

4 

7.i?j6a-oc-i;i j. 

0. 

0* 


0. 


0. 


FIELD angle : J* 
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local oscillator CPTICAu oata« 

N -3 1.390CO:j j£*ac H 0. P^'RO 0, 

M LQ XLINV 0. MAG 0. IFLGl 0 

surface oftfttMiTTERS 

SURF. CURVATURE ScPARATIOM INOEX O3SCURATI0M OATA VIGNETTING DATA 

NO. RwoS T3 THUS BETAP* BETAPY BETAS* BtTASY 

II) c. 0 . l.COOOOJ c. U. Q. 

11 6 . 99 (^a-j 5 d£-^? l.ll^JuJljE♦^^ i.caooai] J. c. C. 0 . 

i,z c. c. L.LOdaaa )• o. o. o, 

special surfaces 

surface RMOl RH02 CRS CC PHl AHP FREQ 

ii L. J. •i.caovOiid j.ooocooo 0. joojQw. o.ajoo.iao o.aiooooo 

(ANGLES INPUT GREATER THAN 6.263 ARE ASSUNEO TO BE IN DEGREES) 

FIELD angle s C.. 

NXY « 25 NOET s 1 mFRS = 26 

OEUFS = 20HEGA = C. VLAHOA = 1. ; 61 1 365 JE-0 2 

detector OATA. 

detector no. 1 I06TEC « 1 J 0 0 0 
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PARAXIAL PAY TRACE 


RECEIVED signal OPTICS 


A 

1.21212121E-32 
1.22025129E-03 
1.2202S129E-33 
-7.0>*723a06E*a2 
-7.0‘»723«06E-02 
-7.04723a06E-i32 
-7.04723906E-32 
-7.0-*723806E-02 
-7.0%723806E-32 

PARAXIAL ENTRANCE PUPIL POSITION T ( 0 ) = » 3 . 7 U 92742E* C 3 TEXIT= -9. S 924922&E* 3 1 

INVERSE OPJECT DISTANCE X.INV = u. 

focal length s -i,32P79b33t*03 9ACK FOCAL LttiGTH = 5,62359C‘^6E + 01 

TS(N-i) = p.fc7750;.CjE*:i Fl = - 1. 320 7956fcE^ )3 


BETA B alpha 

6.25-.. Cw:3E4-C1 3. 0. 

8.2‘>C. uL- 3£ + ':i -2. J1219512E-01 4,46939 594E^ 00 

6.2S.;. JCuuE^JI -2.312X9512E-01 <«,46939394£«^CD 

-l.o69b2747E«:i 6.2HS23HQ2E-02 S, 0 6307 03 3£^C0 

-9.afe67 6525E*wt 6. 2-*6234e(2£- 02 -5,2l6C?6l0E^00 

•3. 4X89285 3E«Cl o. 2 h 623492£- 02 •l«2lS22 6C2E«^0i 

-3.3e7 7-.22 jE*0C 6. 24523492E-02 -1.21875 IC9£*01 

-3.5X2b?f 9«iE«Cl 6. 2<*623492£-02 -1.2046566l£«01 

3.3t729637£-J2 6. 24b23492£-02 - 1.6 3476 355E^01 
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PftRAXIAL pay trace 


LOCAL OSCILLATOR OPTICS 


BETA 

l,390..CLuJE*yt 
l.iSbc JCo.^E«'uC 
l.4ai^664Te-l‘» 


9 

]. 

-1. 25 JOlQCQE-02 
-1.25J0ajCQc«02 


PARAXIAL entrance PUPIL POSITION T(0) = -0. 
INVERSE OPJECT OISTAMCE XLINV * 0. 


ALPHA 

0. 

j • 

9.acooao-.o£«-ci 

TEXiTa 0. 


FOCAL LcN«TM 


l.ll2L^v«b:cYJ2 


SACK FOCAL length s 1* 1120 0 « ) wE^-O 2 


TS<N-1) = 1. 112JCH-E^32 


FL s I*ll2a;00t£»(i2 


A 

2.19424460E-31 

7.19424460E-01 

T.l9424460E-i)l 
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RECEIVER test CASE *1 


CHIEF RAY data for RECEIVED SIGNAL. 

FIELD ANGLE (RADIANS! « j.OOO.O^jO 

CHIEF RAY COORDINATFS ON SJRFACE N 

> Y 7 

0. l.l&u2539tE«Q0 C* 

DIRECTION COSINES ON SURFACE N-l 

L M N 

0. -7.1l94E73oE- 14 l.QICOOOOCE^OO 

RADIUS RW And POSITION TR 3F EXIT PUPIL 

RM TR 

1.536999E3E4^Q2 -9. 6924 3229c^0 1 

TILT AND DISPLACEMENT OF KIT PUPIL IN X- ANO Y-OIRECTIONS 

xtilt ytilt XOISP 

0. -7. 1-542736E-14 C. 


wavefront error iNFORiATION 


VARIANCE s 3.A<*39aC49£- 
RHS s 7.a7i.6d726E-C2 
WAVELENGTH s 1 . 06 11385 JC- J2 
MAXIMUM UNN0PMALI2E0 £RROR> 0. 

MINIMUH UNNORMALI2EO ERRORS- 1 • 64268496E- 0 3 
MAXIMUM normalized ERRORS 
MINIMUM normalized ERRORS • i , S l634 0 HE- 0 1 
APPROXIMATE STRIHL PATIO * 7 , 55 19137 l£- C 1 


to ISP 

1.16025394E^10 
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Reproduced from 
bes( available copy. 



RiCcIV£=> ilCT cas- *l 


aSF {PefiD, FOF RtCiWcO StGNA-.. 

QUflORahT i, ICC, (ORRiS^ONOa TO 1. 3C ‘*<*9t5tE*j3 


L 

1 

1 



11--- 

-1 -1 

-2 

-3 -3 -4 -4 

-5 -5 

-5 

-5 

-4 

-4 

-3 

- 2 

-2 - 1 

2 

2 

1' 



2 2 1 1 

1 u 

-1 

-1 -2 -3 -3 

-4 -4 

-5 

-5 

-5 

-4 

-4 

- 3 

-2 -2 

$ 

7 

r 


3 

3 3 2 2 

2 1 

1 

-1 -1 -2 -2 

-3 


-5 

-5 

-5 

-4 

- 4 

-3 -2 

2 

2 



3 

3 3 3 3 

3 2 

2 

-9 

-2 -3 

-4 

-4 

— — 

- 5 

-5 

- 4 

-4 -3 

1 

1 

1 


2 

2 6 3 3 

3 3 

2 

2 1 1 -u 

-1 -2 

-3 

-3 

-4 

— 4 

-5 

-5 

-4-4 

-1 

-1 

- j 

-c 


1 1 2 > 

3 3 

3 

3 2 2 1 

- -1 

-2 

-3 

-.3 

-4 

-5 

-5 

.4 

-3 

-3 

-3 

-2 

- j 

-1-1 : 1 

2 2 

3 

3 3 2 2 

1 : 

-1 

-2 

-3 

-4 

-4 

- 5 

-5 -5 

•6 

-o 

*6 

- V 

-5 

-4 -3 -2 -1 

•) 1 

2 

2 3 3 3 

2 1 

J 

-1 

-2 

- 3 

-4 

-4 

-5 -5 


-9 

-9 

-ft 

-ft 

-7 -6 -5 -3 

-2 -1 


12 3 3 

3 2 

1 

S 

-1 


-3 

- 4 

-5 -5 

12- 

12- 

11- 

1*- 

1 

-9 -d -7 -3 

-H -I 

-1 

i 1 2 3 

3 3 

2 

1 

j 

-1 

-2 

- 3 

-4 -6 

1**- 

12- 

II- 

1:- 

13- 

12-U-lG -ft 

-7 -5 

-3 

-2 -a i 2 

5 3 

3 

2 

1 

-3 

-1 

- 5 

-4 -5 

lu- 

14- 

!••- 

14- 

14- 

13-13-12-U 

-9 -? 

- ft 

- ,4 - 2 “ j 1 

2 3 

3 

3 

7 

1 

-1 

-2 

- T - w 

n- 

13- 

1^- 

l3- 


14-14-13-13- 

ll-l- 

-6 

-6 -4 -2 G 

2 3 

3 

3 

3 

7 

Cr 

- 1 

-2 -4 

1 j- 

Iti- 

i:- 


12- 

13-13-1^-14- 

13-12- 

Ig 

-ft -6 -4 -1 

: 2 

3 

3 

2 

7 

1 

- 3 

-i -1 



-5 

-f 

-ft- 

10-12-13-14- 

14-13- 

12- 

n -6 -6 -3 

-i i 

2 

3 

3 

3 

2 

1 

-1 -2 


3 

9 


- j 

-5 -6-11-12- 

m-14- 

13- 

12-10 -7 -5 

-2 L 

2 

3 

3 

5 

3 

1 

-i -2 

13 

13 

11 

0 

5 

1 -3 -7-1,- 


1 4- 

13-12 -9 -7 

-4 -1 

1 

3 

3 

4 

3 

2 

: -i 

2? 

24 

£ 2 

lu 

13 

8 3 -2 -7- 

11-13- 

t-4- 

14-13-11 -8 

-5 -2 


2 

3 

4 

3 

2 

1 - . 

37 

3fc 

34 


13 

17 1 J 3-3 

-ft-l2- 

14- 

15-14-12-lG 

-7 -4 

-1 

1 

3 

4 

4 

3 

2 t 

Sj 

49 

46 


34 

26 Id 9 i 

-5-1^- 

13- 

15-15-13-11 

-8 -5 

-2 

1 

3 

4 

4 

3 

2 0 

6i 

62 

5n 

6 2 

>•4 

35 25 10 7 

-t -7- 

12- 

14-15-14-12 

-9 -6 

-3 

0 

? 

3 

4 

3 

2 1 

7«> 

74 

7. 

63 


23 62 It 

2 -5- 

ll- 

1-.-15-15-13 

-1 :• -7 

-3 

- j 

2 

3 

4 

4 

3 1 

ae> 

64 

79 

/ 


51 23 27 i5 

5 -3 

-9- 

13-15-15-14 

-11 -7 

-4 

-1 

2 

3 

4 

4 

.3 1 

93 

92 

87 

79 


57 4, 31 19 

a 

-6- 

13-15-15-14-11 -8 

-5 

- 1 

1 

3 

4 

4 

J 2 

9d 

97 

92 

83 

73 

bg 4? j3 2 1 

3 -c 

-6- 

1 i-l5-lb-l4 

-12 -8 

-5 

-2 

1 

3 

4 

4 

3 2 
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quadrant 2, IJO CORRESPONDS TO 1* 3C«»491SLE4^ J3 

-I -2 -2 -5 -^4 -*» -5 -5 -5 -5 -4 «4 -3 -3 -2 -1 -1 -3 3 1 1 1 1 1 1 

-2 -2 -3 -4 -4 -S -p -S -4 -4 -3 -3 »2 “1 -I -0 1 I 1 2 2 2 2 2 2 

-2 -3 -4 •*, -5 -5 -5 -4 -4 -3 -2 -2 -I -0 1 l 2 2 2 3 3 T 3 3 3 

•3 -4 -4 -5 -S *4 -4 -3 -2 -I -u 0 1 2 2 2 3 3 3 3 T 3 2 2 

-4 -4 -S -4 -4 -4 -3 -2 -1 -u 1 ' 1 2 2 3 3 3 3 2 2 2 2 1 1 

•4 -5 -p -4 -3 -3 -2 -I 0 I 2 2 3 3 3 3 2 2 1 1 3 -0 -0 -1 

-5 *5 -5 -4 -4 -3 -2 -i J 1 2 2 3 i 1 2 2 I 0 -1 -1 -2 -2 -3 -3 

-5 -5 -4 -4 -i -2 -i 0 I 2 3 3 3 2 2 l 3 -1 *2 -3 » -5 -S -6 -6 

-5 -5 -4 -3 -2 -I 0 I 2 3 3 3 2 I 1 -1 -2 -3 -4 -6 -7 -7 -6 -9 -9 

-5 -4 *3 -2 -i 3 1 2 3 3 3 2 1 0 -1 -3 -4 -6 -7 -8 -9-13-11-11-12 

-5 -4 -3 -2 -j I 2 3 3 3 2 I -J -2 -1 -5 -7 -6-10-11-12-12-13-1 3-13 

-4 -3 -2 -i I 2 3 3 3 2 1 -0 -2 -4 -6 -6 -9-11-12-13-13-14-14-14-14 

-4-3-1 i. 1 2 3 3 3 2 u -2 -4 -6 -6-10-11-12-15-14-14-14-13-1 3-13 

-3 -2 -- i 2 3 3 3 2 1 -1 -3 -6 -8- 10- 12- 13-13- 14-13-1 3-12- 11-1 0- 1C 

-2-1 C 2 3 3 i 3 I -1 -3 -5 - 5-1 0- 12-13-14-1 3- 13-12-10 -8 -7 -5 -4 

-2-0 1 2 3 3 3 2 j -2 -S -7-10-12-13-14-13-12-11 -8 -6 -3-0 2 3 

-1 0 2 3 3 3 3 1 -1 -4 -6 -9-11-13-14-14-12-10 -7-3 1 4 811 13 

-1 1 2 3 4 3 2 C -2 -5 -8-11-13-14-14-13-11 -7-3 3 ft IT 18 22 24 

-a 1 3 4 4 3 2 -1 -5 -6 -9-12-14-14-14-12 -8 -3 3 lO 17 23 29 33 36 

0234431-2-5 -8-11-13-14-14-13-10 -5 2 9 17 2b 34 40 46 49 

1 2 34 M 2 J-2-5 -9-12-14-15-14-12 -8 -1 6 15 25 35 44 52 58 62 

1 24 4 3 2 -J -3 -3-10-13-IH-15-14-11 -5 2 11 21 32 43 53 62 o9 74 

134432-1-4 -7-10-13-15-15-13 -9 -3 5 15 26 38 50 62 72 79 84 

134431-1-4 -5-11-14-15-15-13 -9 -2 7 18 30 43 56 68 79 37 92 

134431-1-5 -5-11-14-15-15-13 -8 -1 9 20 53 46 60 72 83 1l 96 
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QUADRANT i, LjO CJRRcS^JNOS TO 1« iC «*4915Ce«^a3 

1 344 3 1-1-5 -i-ll-14-15-15-12 -% -C 9 21 3u 47 61 73 84 98 

1 3 44 3 1-1-5 -i-lt-l4-l5-lS-l3 -9 -1 9 20 33 46 65 72 83 9b 

134431-1-4 -8-11-14-15-15-13 -9 -2 7 18 30 43 56 68 79 37 92 

1 3 - 4 3 2 -1 -4 -<^-lU-t3-l5-l5-l3 -9 -3 5 15 26 38 5'j 62 72 '9 84 

1 2 4 - 3 2 -3-3 -b-lJ-13-14-15-14-11 -5 2 ll 21 32 43 53 62 S9 74 

I 2 3 4 4 2 j -2 -6 -9-12-14-15-14-12 -8 -t 6 15 25 35 44 52 58 62 

] 2 3 - 4 3 1-2-5 -d-ll-13-14-14-13-10 -5 2 9 17 26 34 40 -6 49 

-J 1 54432-1-3-6 -9-12-14-14-14-12 -8 -3 3 10 17 23 29 33 36 

-1 1 4 j 4 3 2 0 -2 -5 -8-11-13-14-14-13-11 -7-3 3 ft 13 18 22 24 

-1 0 2 3 3 3 3 1 -1 -4 -6 -9-11-13-14-14-12-lC -7-3 1 4 8 11 13 

-2 -C 1 4 3 3 3 2 3-2-5 -7-lu- 12- 13-14-1 3-12- 11 -8 -6 -3-0 2 3 

-2-1 - 4 3 3 3 3 I -1 -3 -5 - 9- 1 0- 12- 13-14- 13- 13-12-13 -8 -7 -5 -4 

-3 -2 1 2 3 3 3 2 1 -1 -3 -6 - 8-10-12-13-1 3- 14-13-13-12-11-1 0- 10 

-4-3-1 L 1 2 3 3 3 2 L -2 -4 -6 -9-10-11-12-13-14-14-14-13-1 3- 13 

-4 -3 -2 -I 1 2 3 3 3 2 1 -0 -2 -4 -6 -8 -9-11-12-13-13-14-14-14-14 

-5 -4 -3 -4 -- 1 2 3 3 3 2 1-3-2 -3 -5 -7 -8- 10-11-12-12-13-1 3-13 

r5 -4 -3 -4 -1 u I 2 3 3 3 2 1 0 -1 -3 -** -6 -7 -8 -9-10-11-11-12 

-5 -5 -4 -3 -2 -I 0 1 2 3 3 3 2 1 1 -1 -2 -3 -4 -6 -7 -7 -8 -9 -9 

-5 -5 -4 -3 -2 -1 ^ 1 2 3 3 1 2 2 I u -1 -2 -3 -5 -5 -6 -6 

-5 -5 -5 -4 -4 -3 -2 -1 J I 2 2 3 3 3 2 2 1 0 -1 -1 -2 -2 - 3 -3 

-4 -5 -5 -'V -4 -3 -3 -2 -1 0 1 2 2 3 3 3 3 2 2 1 1 0 -0 -0 -1 

-4 .4 .5 -5 *4 -H -3 -2 -1 -C 1 1 2 2 3 3 3 3 2 2 2 2 I I 

-3 -4 -4 -6 -5 -5 -4 -4 -3 -2 -I -u J 1 2 2 2 3 3 3 3 3 3 2 2 

-2 -3 -H -4 -5 -5 -5 -4 -4 -3 -2 -2 -1 -0 1 1 2 2 2 3 3 3 3 3 3 

-2 -2 -3 -4 -4 -5 -5 -5 -4 -4 -3 -3 -2 -1 -t -0 1 1 I 2 2 2 2 2 2 

-t -2 -2 -3 -4 -4 -5 -5 -5 -5 -4 -4 - 1 - 3 -2 -I -1 -C 0 1 1 1 1 1 1 
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QUADRANT 4, 130 CORRESPONDS TO 1. 3C4491SCE403 


100 

96 

9? 

66 

74 

62 

46 

34 21 lU -7-12-15-16-14-12 

-6 

-5 

-2 

1 

3 

4 

4 

3 2 

96 

97 

92 

83 

73 

60 

47 

33 21 9 -u -6-13-15-16-14-12 

-6 

-5 

-2 

1 

1 

4 

4 

3 

2 

93 

92 

67 

79 

69 

57 

44 

31 19 6 -1 -6-13-15-15-14-11 

-8 

-5 

-1 

1 

T 

4 

4 

3 

2 

66 

64 

79 

72 

62 

5l 

39 

27 15 5 -3 -9-13-15-15-14-11 

-7 

-4 

•1 

2 

3 

4 

4 

3 

1 

75 

74 

7C 

63 

54 

44 

33 

22 11 2 -5-11-14-15-15-13- 

-10 

-7 

-3 

-0 

2 

3 

4 

4 

3 

1 

63 

62 

S6 

52 

44 

35 

25 

16 7 -1 -7-12-14-15-14-12 

-9 

-6 

-3 

0 

2 

3 

4 

3 

2 

1 

SO 

49 

46 

41 

34 

26 

16 

9 2 -5-10-13-15-15-13-11 

-8 

-5 

-2 

1 

3 

4 

4 

3 

2 

L; 
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INTERVAL = 1. u362y5«#5E-J2 NUNBER INTERVALS * 5l SUN = 0. 
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PH4SE 


QUAORSNT I, 10 J COPRESPONJS TO THO PI 

-16-16-l6-ifc-l6-l7-l?-19-J2 tl 39 39 33 38 39 39 39 19 I.J Al 41 43 45 »8-45-35 

-l6-l6-16-lo-l6-16-lb-lb-l3-l7-23 •»S 4u 39 39 39 39 39 40 40 41 42 43 45 48-45 

-17-l7-17-17-17-16-lb-lb-l=)-iei-16-l»-2b 44 40 39 39 19 39 40 40 4l 42 43 45 48 

-19-19-13-la-18-17-l7-l7-li-l6-lb-lS-l6-l9-42 42 40 39 39 39 40 40 41 42 43 45 

-21-21-21-2i.-2J-l9-18-18-l7-l7-lfj-lb-16-16-17-26 44 40 40 39 40 40 40 4l 42 43 

-26-26-25-24-23-22-21-20-10-18-17-17-14-16-16-17-21 48 41 40 40 40 40 -0 41 42 

-36-35-34-33-30-28-26-23-22-20-19-18-17-16-16-16-16-19-48 41 40 40 40 -0 41 41 

-47-47-46-45-42-39-35-31-27-24-21-19-18-17-16-16-15-16-19-48 41 40 40 -0 40 41 

45 45 46 4? 48 50-47-43-37-31-26-23-20-18-17-16-15-15-16-19 48 41 40 40 40 40 

42 42 42 43 43 44 46 4a--.8-43-36-29-24-20-13-17-l6-lS-l5-15-20 45 41 <0 40 40 

40 40 40 4L 41 42 43 44 45 48-46-30-31-25-21-18-17-15-15-15-16-26 43 4 1 40 40 

38 38 38 39 39 40 41 42 -3 44 47-48-41-32-25-20-18-16-15-15-15-16-43 42 41 40 

34 36 36 37 37 38 39 40 -1 42 44 46 50-42-32-24-19-17-15-15-14-15-19 46 41 41 

33 34 34 35 36 36 37 38 40 41 42 -4 46 50-42-30-22-18-16-15-14-14-15-30 43 41 

29 29 3o 3l 33 34 36 37 38 39 41 42 44 46-48-39-27-20-17-15-14-14-14-17 49 42 

21 21 23 25 23 31 33 35 37 38 40 4l 42 44 47-47-35-24-18-16-15-14-14-15-28 44 

9 9 11 14 19 24 29 33 35 37 39 4j 41 43 45 48-43-29-20-17-15-14-14-14-18 48 

1 2 2 4 7 13 21 28 33 36 38 39 41 42 43 46-48-37-24-18-15-14-13-14-15-38 

-1 -1 - 1-4 1 3 9 19 29 34 37 39 40 41 43 44 48-43-28-19-16-14-13-13-14-23 

-3 -3 -3 -2 -I -0 2 8 22 32 36 38 40 41 42 43 40-48-34-21-15-14-13-13-14-18 

-4 -4 -4 -3 -2 -1 2 11 27 35 38 39 40 42 43 45 49-40-23-17-15-13-13-13-16 

-5 -5 -5 -4 -4 -4 -2 -0 4 20 33 37 39 40 41 42 44 47-44-26-13-15-13-13-13-15 

-5 -5 -5 -5 -5 -4 -3 -2 3 12 31 35 39 40 41 42 44 46-46-29-18-15-14-13-13-14 

-6 -6 -b -5 -5 -5 -4 -3 -0 7 28 56 39 40 41 42 43 45-48-52-19-15-14-1 5-13-14 

-6 -6 -6 -6 -5 -5 -4 -3 -1 4 26 36 38 40 41 42 43 45-49-34-20-15-14-13-13-14 
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QUADRANT 

4, 

lu 
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; Tc 

1 TWO 
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42 43 46-48-37-24-18-15-14-13-14-15-38 
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29 

33 

35 

37 

39 

40 
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43 45 48-43-29-20-17-15-14-14-14-18 46 
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25 

28 

31 

33 

35 

3/ 

38 

4C 

4l 

42 

44 47-47-35-24-18-16-15-14-14-15-28 44 

29 

29 

3w 

3i 

23 

34 

3o 

37 

33 

39 

41 

42 

44 

46-48-39-27-20-17-15-14-14-14-17 49 42 

33 

34 

34 

35 

3b 

3b 

27 

38 

4 j 

41 

42 

44 

46 

50-42-30-22-18-16-15-14-14-15-30 43 41 

. 3b 

36 

36 

37 

37 

36 

39 

40 

41 

42 

44 

40 

50- 

42-32-24-19-17-15-15-14-15-19 46 41 41 

38 

36 

36 

39 

39 

4ii 

41 

42 

43 

44 

47- 

48' 

-41- 

32-25-20-18-16-15-15-15-16-43 *♦ 2 41 40 

40 

40 

4u 

41 

41 

42 

43 

44 

4 3 

48- 

46- 

39 

-31- 

25-21-18-17-15-15-15-16-26 43 4l 40 40 

42 

42 

42 

43 

43 

44 

46 

48- 

48- 

43- 

36- 

29 

-24- 

20-l8-l7-16-l5-l5-l5-2o 45 41 >4 0 40 40 

45 

45 

46 

47 

48 

50- 

47- 

43- 

37- 

31- 

26- 

23 

-20- 

18-17-16-15-15-16-19 48 41 40 40 40 40 

-47- 

47-46- 

4P- 

42 

-39- 

35- 

31- 

27- 

24- 

21- 

19 

-18- 

17-16-16-15-16-19-48 4i 40 40 40 40 41 


-3b-35-3‘.-3J-20-28-i6-23-22-20-l9-ia-17-16-16-16-16-l9-<»8 41 4j 4) 4C 40 41 41 

-26-26-25-24-2J-22-21-2J-19-18-17-17-16-16-16-17-21 48 41 40 4) 40 40 40 41 42 

-2l-21“2l*2w“20-lS-l8-l8-l7-l7-lfe-lb“ib*i6*l7-26 44 40 40 39 4o 40 40 •» 1 42 45 

-19-l9-lo-ib-l8»l7-l7-l7-l3-l6-l6-l6-16-l9-42 42 40 39 39 39 4: 40 41 ^2 43 45 

-l7-i7-17-l7»l7-16-l6-lb»lD-l6-lo-l8-26 44 40 39 39 39 39 40 4C 4l 42 ^3 45 48 

-lb-l6-lb*lc-lb-16-lb-16*l3*l7-23 <»5 40 39 39 39 39 39 40 40 41 42 43 w5 48*45 

-lb-lb-lfe-lb-lb-l7-l7-l9-32 43 39 39 38 38 39 39 39 39 40 41 41 43 45 48-45-35 

FOR THIS Flttr 4NGLE ANO WAVELENGTH 26 RAYS WERE VIGNETTED OR 09SCUREO, 

c rats failed to converge in SSRT. 




154 


RECEIVER TEST CASE *1 


CHIEF rat data for local DSCILLATOR. 


FIELD ancle (RADIANS) ^ U.QuaJOJJG 
CHIEF ray coordinates ON SJRFACE N 
X Y 

0. a. 

direction cosines on surface n-1 

L i 

0. u. 


0. 


z 


N 

i.ooooooacE>oo 


RADIUS RH and position TR 3F EXIT P'JPIL 

RH TR 

i.ii20(^:gqe«'OZ u. 

TILT AND OISPLACEMENT OF EXIT PUPIL IN X- AND Y-OIRECTIONS 

XTILT YTILT XOISP 

0* i. 0. 


1 . 


YD ISP 


MAVEFRONT error INFOR.iATION 


VARIANCE « 3.6331S7Q6E-1S 

RNS a L.7lbu9a76£-16 
wavelength 3 I«a61l3050£-J2 

HAXIKUH UNNORHALIZEO ERROR* S • 59 1 USIZE* l7 
HINIHUH UNNORMALIZEO £RR0R*-2« Z2Qv4605E-l6 
HAXIHUN NORNALIZEO ERROR* 5 • 23l282i7E*tS 
HINIHUH NORMALIZED ERR0R#«2. C 925 1295E* IL 
APPROXIMATE STREHL RATIO * I « 00 900 J C0£4>0 0 
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RECEIVER TEST CASE *1 


ASF (REAL). FOP LOCAL OSCILLATOR. 

QUADRANT 1. luO CORRESPONDS TO a. UC4360Q3e«QG 
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QUAORANT 2, IJO CORRESPONDS TO 2. U( 4 36003C^UQ 
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quadrant 3, 100 CORRESPONDS TO 


2« Jb436ad3£4d0 
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69 
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quadrant 4, lOQ CORRESPONDS TO 2. >lI.43b()a3E«aC 
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RECEIVE® TEST CASE *1 


ASF <IMACINAPV» ,FOF LOCAL OSCILLATOR. 

OUAORArJT 1, 1U3 CORRESPONDS TO 2.H204e796£*09 

-7l-7l-7:-7t-69-b6-€ 7-65-bi-6l-57-53«5J-45-42-3a-35-30-22-ie-14 -9 -5 -3 -I 0 
-7l-7l-7l-7u-7 j-69-b7-bt»-64-bl-57-53«49-45-41-36-34-29-2l-l7-13 -« -4 -2 *0 1 

-7l-7l-7l-7t-7D-b9-67-b5-*>3-bl-57-53-.49-44-4l-37-33-2«-23-lb-U -6 -2-0 2 3 

-70-70-7 j-7v-b9-b0-b7-b5-63-bO-5b-52-48-43-39- 36-32-27-19-14 -9-4-0 2 4 5 

-67-b7-bb-6o-b>-6».-63-61-39-5b-52-48-4-*-39-3b-32-2«-23-15-10 -50 4 6 8 9 

-6l-bl-bl-b--e.J-59-57-55-53-51-47-43-39-34-30-27-23-l7 -9 -5 -0 5 9 11 13 14 

-57-57-56-56-65-54-53-51-^9-46-42-36-34-29-26-22-18-13 -5 0 5 10 14 16 17 16 

-52-52-52-52-61-50-49-47-43-42-38-34-30-25-21-16-14 -6-0 $ 9 15 19 20 21 22 

-43-43-43-42-42-41-40-38-36-33-29-25-21-16-13 -9 -5 -0 8 13 17 23 27 28 29 30 

-37-37-37-3/ -36-36-34- 32-3u-26-24-2u-16-11 -7-4 0 5 14 18 23 28 32 33 34 35 

-34-34-34-33-33-32-21-29-27-24-20-16-12 -7 -3-0 4 9 18 22 27 32 36 37 38 38 

-30-30-3.-31-29-29-27-25-23-21-17-12 -6 -4-0 3 7 13 21 26 30 36 39 4l 41 42 

-27-27-27-2C-26-25-24-22-2 j-17-13 -9-5 0 4 7 11 16 25 29 34 39 43 44 45 45 

-22-22-22-21-21-20-19-17-15-12 -8 -4-0 5 8 12 16 21 30 54 39 44 46 49 49 50 

-13-18-16-1^-17-16-15-13-11 -8-4 0 4 9 12 16 20 25 34 38 43 48 52 53 53 53 

-14-14-13-13- 13-12-11 -9-7-4 C 4 8 13 17 20 24 29 36 42 47 52 56 57 57 57 

-9 -9 -9 -9 -8 -8 -6 -4 -3-0 •* 8 12 l7 21 24 28 33 42 46 51 56 60 5l 61 61 

-7 -7 -6 -c -6 -5 -4 -2 -J 3 7 11 15 20 21 27 30 16 45 49 53 59 63 63 64 63 

-5 -5 -5 -4 -4 -3-2-0 2 4 9 13 17 22 25 29 32 38 47 5l 55 61 65 o5 66 65 

-3 -3 -3 -2 -2 -1 -J 2 « 6 11 15 19 24 27 31 34 40 49 53 57 63 67 67 67 67 

-2 -2 -1 -I -1 -li I 3 5 8 12 16 20 25 29 32 36 41 50 54 59 64 68 69 69 68 

-I -1 -I -g 1 2 4 a 8 13 17 21 26 29 33 36 42 51 55 60 65 69 70 70 69 

-J -C -: -g J 1 2 4 6 9 13 17 21 26 33 33 37 h2 52 56 60 66 70 70 7 j 70 

-0-0-0 L 1 1 3 5 5 9 13 19 22 27 30 34 37 43 52 56 61 66 70 7i 7l 7£ 

-0-0 u g 1 2 3 6 7 9 14 18 22 27 30 34 37 «*3 52 57 61 67 70 7l 71 71 
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QUADRANT 2f lOd CORRESPONDS TO 2. 42 0487<3flE*a9 


-e9-9l-93*95.-9**-92-91*90-83-a3-52«d2-8a*7 8-77-75-73-r2-72-/l-7l-7 0-71-7 i-7l 
.9l.94.9S.97>9b-94-93-92-d/«8$-a4-9 3-81-79-78-7b-7<**73-72-7l-7l-7l-7l-7 1-7 1 
-93.95.97-99-98-96-94-93-99-86-84-83-82-80-78-76-74-73-72-71-71-71-71-71-71 
-95.97-99*0 6- 99-97-95-94-89-86-95-84-82-80-78-76-73-72-72-71-70-70-70-70-70 
-94-95-98-99-98-96-94-92-87-84-92-91-79-77-75-73-70-69-68-67-67-67-66-97-67 
-92-94-90-97-96-93-91-69-83-80-78-77-75-72-70-68-65-64-63-62-61-61-61-6 1-61 
-91-93-94-96-94-91-89-67-80-77-75-74-72-69-66-64-61-6C-59-58-57 -57-57-57-57 
•9j.92.93. 94-92-89-67-84-7 8-7 4-72-7 1-68-65-6 3-60-57-55-54-53-53-52-52-5 2-52 
-36-87-88-89-87-8 3-6 0-78-7 1-67-65-63-6 1-57-54-52-48-47-46-44-44-4 3-43- •♦3-43 
-83-85-86-86- 04-80- 77-74-67-6 3-01-59-56-53-50-47-43-42-40-39-38-38-38-3 7- 17 
-82-84-84-85-62-78-75-72-55-61-58-56-53-50-47-43-40-38-37-35-35-34-34-54-34 
-82-d3-83-8‘.- 91-77-74-70-63-59-56-54-51-47-44-41-37-35-34-32-31-31-31- 5 0-30 
-80-82-82-82-79-75-72-68-61-56-53-51-48-44-41-37-34-32-30-29-28-27-27-27-27 
-78. 79. 80 -8U-77-72-69-65-57-53-5C-47-44-4 0-36- 13-29-27-26-24-23-22-22-2 2-22 
-77-78-78-78-75-70-66-63-5*-50-47-44-41-37-33-29-25-23-22-20-19-l8-18-l8-l8 
-75-76-76-76-73-69-64-6fl-52-47-^3-«*l-57- 33-29- 26-22-19-19-16-15-14-14-1 4-14 
-73-74-74-73-70-65-61-57-48-43-^0-37-34-29-25-22-17-15-13-12-11-10 -9 -9 -9 
-72-73-73-72-69-64-60-55-47-42-38-35-32-27-23-19-15-13-11 -9 -8 -7 -7 -7 -7 
-72-72-72-72-68-63-59-54-46-40-37-34-30-26-22-18-13-11 -9 -7 -6 -6 -5 -5 -5 
-7l-7l-7l-7l-b7-b2-58-53--*4-39-35-32-29-24.2Q-16-l2 -9 -7 -6 -*♦ -4 -3 -3 -3 
-71-71-71-7^-67-61-57-53-44-38-35-31-28-23-19-15-11 -8 -6 -4 -3 -2 -2 -2 -2 
-71-71-71-71-67-61-57-52-4 J-38-34-SI-27-22-18-14-10 -7 -6 -4 -2 -2 -I -1 -1 
-71.71.71-7l-66-61-57-52-*3-38-34-31-27-22-18-14 -9 -7 -5 -3 -2 -1 -I -0 -0 
-71-71-71-70-67-61-57-52-43-37-34-30-27-22-18-14 -9 -7 -5 -3 -2 -1 -0 -0 -0 
-71-71-71-70-67-61-57-52-43-37-34-30-27-22-18-14 -9 -7 -5 -3 -2 -1 -0 -0 -0 
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QUADRANT 3* 100 CORReS>’3NDS TO 2. <*2046798£-09 


-71-71-71-7L-67-61-57-52-U1-37-J4-30-27-22-1S-14 -9 -7 -5 -3 -2 -1 -0 -0 -0 

-7l-7l-7l-7w“fc7-bl-57-52-H3*37-34-3ij-27-22-t8-14 -9 -7 -5 -3 -2 -I -0 -0 -0 

-70-71-7l-7C-fc6-61-56-52-43-37-Jv-30-27-22-l8-13 -9 -6 -S -3 •! -1 -0 -0 0 

-70-70-7D-7L-fcfa-60-56-52«42'37-33-30-26-21-17-l3 *9 -6 •% -2-1-0 0 0 0 

-69-70-7C-69-65-60-55-5l-**2-36-33-29-26-ei-l7-13 -3 -6 -4 -2 -l 0 0 1 1 

-6e-69-69-6e-6‘*-59-55-5u-4l-36-32-29-25-20-l6-l2 -8 -5 -3 -1 a 1 1 1 2 

-87-67-67-67-fc3-57-53-49-4U-34-3l-27-2^-19-15-ll -6-4-2012233 
-65-66-65-65»-fcl-55-61-47-33-32-29-25-22-l7-13 -9 -4 -2 0 2 5 4 4 5 5 

-63-64-63-63-59-53-49-45-36-30-27-23-20-15-11 -7-302456667 
-61-61-61-60-56-51-46-42-33-28-24-21-17-12 -8-4034688999 
-57-57-57-66-52-47-42-38-29-24-20-17-13 -8-4-0 4 7 9 11 12 13 13 13 14 

-53-53-53-52-48-^3-28-34-23-20-16-12 -9 -4-0 4 8 11 13 15 16 17 17 18 18 

-50-49-49-46-44-39-34-30-21-16-12 -6 -5-0 4 6 12 15 17 19 20 2l 21 22 22 

-45-4S-4‘.-H3-39-34-29-25-l6-ll -7 -4 -J 5 9 13 17 20 22 24 25 26 26 27 27 

-42-41-41-39-36-30-26-21-13 -7 -3 0 4 8 12 17 21 23 25 27 29 29 30 30 JC 

-38-38-37-36-32-27-22-18 -9 -4 0 3 7 12 16 20 24 27 29 3l 32 33 33 34 34 

-35-34-33-32-28-23-18-14 -5-0 4 7 11 16 20 24 28 30 32 34 36 36 37 37 37 

-30-29-28-27-23-17-13 -8 3 5 9 13 16 21 25 29 33 36 38 40 41 42 42 43 43 

-22-21-20-19-15 -9-5 0 8 14 16 2l 25 30 34 36 42 45 47 49 50 5l 52 52 52 

-18-17-16-14-10 -5 -0 5 13 18 22 26 29 34 38 42 46 49 5l 53 54 55 56 56 57 

-14-13-11 -9-5 0 3 9 17 23 27 30 34 39 43 47 51 53 55 57 59 60 60 61 61 

-9 -6 -6 -4 -0 5 10 15 23 28 32 36 39 44 48 52 56 59 61 63 64 65 66 66 67 

-5-4-2 L •• 9 14 19 27 32 36 39 43 48 52 56 60 63 65 67 68 69 70 7 0 70 

-3-2 u 2 6 11 16 20 28 33 37 41 44 49 53 57 61 63 65 67 69 70 70 7 i 71 

-1 0 2 4 8 13 17 21 29 34 36 4l 45 49 53 57 61 64 65 67 69 70 70 7l 71 

-0 1 3 5 9 14 18 22 30 35 38 42 45 50 53 57 61 63 65 67 68 69 70 7 0 71 
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quadrant iuO CORPES=^3NOS TO 2. 42S487Q6E-09 
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79 
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27 

28 

29 

30 
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34 

17 
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44 

48 

51 

53 

56 

61 

68 

72 

75 
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82 

82 

8c 

30 

3(i 

3c 

3: 

31 

31 

32 

34 

35 

37 

4l 

44 

47 

51 

54 

56 

59 

63 

73 

74 

77 

81 

84 

8 3 

83 

82 

34 
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34 
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3h 

35 

35 

37 

38 

40 

43 

47 

5 3 

53 

56 

56 

61 

65 

72 

75 

78 
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85 

54 

84 

82 

37 

37 

37 

3b 

3fl 

38 

39 

4u 

42 

43 

47 

50 

5 3 

56 

59 

61 

63 

67 

74 

77 

8j 

84 

86 

56 

85 

63 

43 

43 

43 

4i 

43 

44 

44 

46 

47 

4d 

52 

54 

57 

61 

63 

65 

67 

71 

7ft 

80 

63 

87 

89 

5 8 

87 

86 

92 

52 

52 

5; 

52 

53 

53 

54 

55 

57 

60 

63 

65 

68 

71 

72 

74 

78 

84 

86 

69 

92 

94 

9 3 

92 

90 

57 

57 

57 

5? 

5/ 

57 

16 

59 

6 J 

61 

34 

56 

69 

72 

74 

75 

77 

80 

87 

99 

91 

94 

95 

94 

93 

91 

61 

61 

61 

61 

61 

o 1 

62 

63 

34 

65 

3d 

70 

72 

75 

77 

78 

8J 

83 

89 

91 

93 

96 

97 

96 

94 

92 

67 

67 

67 

6b 

67 

67 

67 

6d 

69 

70 

73 

75 

77 

79 

81 

82 

64 

87 

92 

94 

96 

98 

99 

98 

96 

94 

70 

7u 

7c 

7. 

7w 

70. 

71 

72 

72 

73 

76 

78 

80 

82 

84 

85 

86 

89 

94 

95 

97 

99100 

99 

97 
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71 

71 

71 

71 

71 
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71 

72 

73 

74 

76 

78 

63 

82 

83 

84 

86 

88 

93 

94 

96 

98 

99 

97 

95 

93 

71 

71 

71 

7i 

71 

71 

71 

72 

7 3 

74 

76 

78 

79 

61 

83 
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85 

87 

92 

93 
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94 

91 
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71 
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71 
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72 

72 

73 

75 

77 

78 

80 

82 

82 

83 

86 

90 

91 

92 

94 

95 

93 

91 

69 

INTERVAL 

= 

1. 

. :3629545i 

>32 

NUMSEr^ 

[NTERVAL5 

a 51 SUM 

= 

0. 
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phase 


QUADRANT 1, ijC CORRESPONDS TO TMO PI 
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QUADRANT 4, lOJ CORRESPONDS TO TMQ PI 
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RECEIVER TcST CASE *1 


DETECTOR CO^’PUT AT IONS . 


DETECTOR NO. 5 1 

RECEIVED SIGNAL BEAM IS SMIFTEO X = 0. AND V * ii , 

WITH RESPECT TO CEfTER 0? DETECTOR. 

LOCAL oscillator BEAM IS SHIFTED X = C. AND t = 0. 

WITH respect rO CENTER OF DETECTOR. 

CIRCULAR OLTECTOR DIAMETER = 2. C T26qd9lE- 01 

SIGNAL ®OHtO = l.e947bo3I^-<* 

LOCAL OSCILLATOR POWER - I . 285 i427E-0 1 
CROSS PRODUCT POWER = 3. a J 2? A26E-f J1 

PHASE MATCH zFFICIEfCY = 9.S95U11E-01 ( -.018 09 I 

average phase SHIFT (RADIANS) > 2. b263201E-01 

OPTICAL TRAnSMISSIO)* s 9. 7581079E-01 < -.106 09 ) 

FOCUSING EFFICIENCY s 7 . b 43 1 7 2b£- 1 1 ( -1.167 OB ) 

NON-HETERCOtHt DETECTION EFFICIENCY s 7 . 45 629u3E-31 ( -1.274 OB 

L.o. Illumination efficiency a 6,29ii«354£-oi ( -.ai4 OB ) 

MAXIMUM AhTtNNA GAIN a 2 . 38 o29 15E 9 9 ( 93.777 OB ) 

RECEIVER fcFFICIENCY TO I.F. * 4. 72J4378E-L1 i -3.2b0 OB ) 
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APPENDIX D 


■^JCEIVE^. Ti'^T CAS£ *2 

SIMP BITAtdl = sJ.irS, 2.621, M(ll = t.,.., IPLGltl) = 0,0, 

N(l> = Mdl = I.IJ, tLINlMlI = 2.. -9. *-35695 7-4, 

9A0Sd) = -1026., J., 52.61152, 3-,. 0866, 3*3.. -105.120**, -161.596*, 

940Sdl) = 3«„., 1163.18356, 666. 2130-., 0., -105.120*-, -161.569*, 0., 3., 
rsd) = C., -51*., -32.1, -2., -910., t., 0.. -193., -2..0., 

TSdll = -160.3, 0., u., -2., -392., -193.. -2.,i)., -100.3, J.. 

XiUSdl =_ 1 . , -1., -1., -,.30062, *•-!., -*.10062, ,*-l., 

XMIlSd*l=-*.i -i 162, -1., -1., -*,01062, 3»-l., 

S<aP6<2) = ... 

PZEROdI = -1.-6, 1.-3, 

20MEGA = 1.-J25, 

s-:mo 
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LACOHfl - LASER C3MHUNICAT0R ANALYSIS PROGRAM 


PAGOS CORPORATION 


RECEIVER TEST CASE *Z 


RECEIVED SIGNAL OPTICAL DATA. 


N 

-12 

BETAO 

5.39750MjdE4-0i 

H 0. 

P7ER0 -i.Ou OOOOOCE-06 

H 

1 

XLINV 

0. 

HAG 0. 

IFLGl 0 


SURFACE 

parameters 







SURF. 

CURVATURE SEPARATION 

INDEX 


OBSCURATION DATA 


VIGNE 

TTINC DATA 

NO. 

rhos ts 

XMUS 


BETAPX BETAPY 


BETAS* 

BETASY 

1 

0. 0. 

l.COQOOO 

0* 

0. 

0. 


0. 

2 - 

9.727&26SE- JA-5.1AJj0u0E>02 

-l.COOdOO 

0. 

0. 

G. 


0. 

3 

0. -3.2lj0J00E>01 

-i.LOOCOQ 

3. 

0. 

0. 


0. 

A 

1.693«64 9E-02-2.uu3J0uQE«>i]0 

-A.cao62a 

0. 

0. 

c. 


0. 

S 

2.933b669E-y?-9.U«:JQb0E402 

-i.COOOOO 

u. 

c. 

0. 


0. 

6 

0. 0. 

•i« ooocoo 

0. 

0. 

0. 


0. 

7 

0. 0. 

-l.COOQOO 

0. 

0. 

c. 


0. 

6 

0. -1.9Jw0uQ0E^02 

-1. OQQOOa 

0. 

0. 

0. 


0. 

9 - 

9.Sl26 976E-C3*2.0iiyJuOOE«-QO 

-4.C00620 

0. 

0. 

u . 


0. 

10 - 

b. Idb2967c-J3 0. 

-I.COOOOO 

0. 

0. 

0. 


0. 

11 

0. -1.0034u00E^02 

-I.COOOOO 

0. 

0. 

0. 


0. 

12 

0. Q« 

-i.COOOOO 

0. 

c. 

0» 


0. 

















SURFACE 

SKAPA ALPHA 


BETA 

GAMMA 


BELT# 

EPSItON 

2 

C. 0. 

0. 


0. 

0» 


0. 


field angle s 
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LOCAL OSCILLATOR CPTICA-. OATA 


N -8 

9LTAj 2,i21C3. j JEfQL 

H Q. 


PZERO 1.0. 

ocooo 

u£-03 


rt X3 

XLINV -9#4j5695^ jE-U«* 

MAG 0. 


IFLGl 0 




SURFACE 

PARAMETERS 







SURF, 

CURVATURE SEPARATION 

INQEX 

OBSCURATION DATA 


VIGNETTING 

OATA 

NO. 

RHOS T5 

XHUS 

9ETAPX 

BiTAPV 


BETASX 

BETASY 

13 

0. 0. 

-l.CiOOOQO 

'3 a 

0. 

j. 


0. 

14 

9.405695 7c-34-2. jC J J*j-JOE«-(iO 

-t».C<)Q62Q 

0. 

0. 

2 , 


0. 

15 

1. 496523 7S-C3-3. 92: JE*>j2 

-l.COOuOO 

3. 

c. 

0 . 


0. 

16 

0, -i.93:ju:oe-»02 

-1.000000 

0. 

0. 

b • 


0. 

17 - 

9.5l<.b97 8E-OZ-2.0ajg30Q£»^0 

-4.CQ0620 

0. 

G. 

c. 


0. 

la - 

6.1892908E-33 0. 

-l.C0Qu03 

0. 

0. 

0. 


0. 

19 

0. * 1 a Q l 3 j w u OE^u 2 

-l.COOOQO 

0. 

0. 

0. 


0. 

2Q 

C. Q, 

-laCaOCQQ 

0. 

0. 

Q. 


0. 


FIELD angle = 3. 

NXT a 25 NOET * 1 NFRS = 26 

OELFS a' Q. ZOMEGA a !• 62530 OOQC«-QO VLAHOA a 6ll3d5CE- 32 

DETECTOR OATA. 

detector no. L tOETEC « 1 3 O 0 0 
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PSRaxlAL RAY TRACE - RECEIVED SIGNAL OPTICS. 


P::TA 

9 

ALPHA 

A 

‘j. 7975 »:l0uE*l?1 

J. 

C. 

1.65270959E-02 

f'. 397i>LL--aE*0l 

-l.05J«9728E-01 

0« 

1.85270959E-02 

2.27377fa75£-13 

-1.070C9726E-01 

9.52292728£^00 

I.d5270959£-J2 

•3, 37«.6l226c*ut 

-2.9656721fiE-ai 

l«Qll7o47tE^01 

5*93W9566CE-01 

-3.5l9w72d>i£» jC 

1. 32l253<«6E-02 

l.d41A34e9E«Ll 

-3.23266946E-01 

B.StWi 5363E^ul 

l.32l253**b£-02 

-2.6375ft572£*D2 

•$.2326694bE-31 

8,5l^33363c»uL 

1.321253C6E.32 

-2*«375fl572E^02 

-3.2326694&E»01 

6«Si.t«3i3D3E^jL 

1. 32l293<.bE*02 

•2.83758572E^C2 

-3.23266946E-01 

U-5W3528£t-l 

•3.02329«49c-0i 

-3.46149093£*02 

9.55741745£*30 

1..9»v32ll5E»Cl 

-9.98720J22E-02 

-3.41371125E^02 

3,21663129E^00 

l,.9C3?ll5t*Jl 

-9.9472-‘]22£-02 

-3.41371125E+02 

3.21863129E^30 

8«dfc0^97u3fjJ 

-9, 9d72Ju22E-02 

•1.8S4240&aE«01 

3.21863129E4-00 

oa-RaxiAL enTf^anc* 

POSITION T(3> = -0. 

TEXlTs -1, 

6u 6 w 9b u E ^ 3 2 


INVERSE OPJtRT distance YLINV = D. 

FOCAL LENGIN = S.l J4,ll3.£*iJZ RACK FOCAL. LENGTH = -1.091718531*02 

TS(N-l) = -i. 091719£3E*d2 Fl = 5. «. C AI.175AE* J2 
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PARAXIAL RAV TRACE 


local oscillator optics 


PcTA 

3 *& 21 (.Uwuii£^JC 

JC 

3 . 62 Ub 9 ej 4 C 4 'UC 
9»u43i.iJE4 8c*)Jl 
5.712J5GJ7E-03 
b.HlI>»i735£-J2 
5. 41141755c- 
-4,46b49255£-ai 


B 

}. 

7. 39722o95£-03 
-<*. 3 a 93183 lE-a 3 
-4.399J1631E-03 
-6. Cl95ai29£-Q3 
-5.<Uh59063£-03 
-5,ait»59053£-03 
-5.0l459<]53£-03 


ALPHA 

0. 3 
J. 3 
1.90737 319E-01 3 
1.49416222E^02 3 
2.22866933E^02 -5 
2.19d96627E*-02 -I 
2.19a96627E^02 -1 
2.95590392E«>01 -I 


PARAXIAL entrance PUPIL POSITION T(0I s -0. 


TEXITs -I.i5d763a5£^02 


INVERSE OBJECT DISTANCE XLINV a -9. 405695 70E-C4 


FOCAL LENGTH a 5. 226747b }E»02 BACK FOCAL LENGTH s -1. 0791 344dE^0 1 


TS(N-l) a -1. 0d5d9d93E4>u2 FL a .3«616643a2E^02 


A 

ai533766E-0l 
81533766E-01 
80677257E-01 
8J677257E-01 
93153915E^J0 
89766283E-*-0 0 
89768283E«^00 
89768283E«^00 


171 



RECtIVEK CASE *2 


CHIEF RAY data for RECtiVEO SIGNAL. 


FIELD AMGLL (RADIANS) s w.uuaLUJOC 
CHIEF ray CDOROINATES ON SJRFAC; N 

X Y Z 

0. J. c. 

OI-tECTION COSINES ON SURFACE N-1 
L H 

0. j. 

RADIUS RH AnO POSITION TR OP EXIT PUPTL 

RW TR 

-3* lUd9Z5SE«L j -1. ub 06096 JE4 02 

TILT ANO OIS®LAC£HENT OF £<IT PUPIL IN X- ANO Y-OIRECTIONS 


N 

-l^OOOOOOCOE^OO 


XTILT 

0 . 


YTILT 


XOISP 

6 * 


YDISP 


•wavefront error information 


variance s 2.3169S636E-S2 
RHS s 3.65t397uuE-£2 
wavelength s I.u6lli650c-J2 
MAXIHUH UNNORHALIZEC ERROR* 0« 

HINIHUM UNNOPNALIZEO ERROR**tt« ll9593<*l£-04 
HAXIMUH NORmALIZEO ERROR* J. 

HINIHUH NORMALIZEO E RRORs* 7 • 69177534E* 32 
approximate STRcHL AATIO s 9. 5V0fll996£-ri 
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OI.ReCTION ON SURFA^£ N-i 


L H N 

0 . 0 . * 1 . 0 j 0 j 0000£«>00 

RflOlUS RH ANO POSITION TR OF EXIT PJPIL 

RW TR 

-r.27l3‘j.j7lE»01 -l.l->a7o385£*Q2 

TIlT ANO OISPLACEMENT op exit PUPIL IN X- ANO Y-OIRECTIONS 

XTKT YTILT XOISP YO ISP 

0* j. 0. o* 


wavefront error inforiation 


VARIANCE s 2.<*2966&7i*£-w? 

RMS = 3.6'*l5‘«9‘*l£-C5 

wavelength : 1 . 0611 J8St}E« 

MAXIMUM UKNOPHALIZEO ERROR* 3 . 72 1N2MN4E* I 0 
MINIMUM UNNOPMALlZSr ERROR** tt . 6 59659C2E* 3 7 
MAXIMUM NCRMAlIZEO error* 3 . 5 Q 7 J U 0 5E* 06 
MINIMUM NORMALIZEO E RROR* - 9 . I6u 72^566- 05 
APPROXIMATE STREhL PATIO * 9 . 999999U6E-0 I 
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receive-? TtST CASE *2 


detector CUMPUTATIONS. 


detector mo. = 1 

RECEIVED liir.NAL IS SHIFTED X s 0. ANO Vs D, 

with respect TO CEMER OF DETECTOR. 

LOCAL OSCILLATOR 9fcAH IS SHIFTED X = C. AND Y = C. 

with respect to center OF DETECTOR. 

CIRCULAR DETECTOR DIAMETER = 1. 29b24fll4£-31 

signal power = 6.21099u5£-l7 

local oscillator power = 1.263 J695E-J5 

cross product power = 2.d425l64£-i)3 

PHASE MATCH EFFICIENCY s i. 962 135 5E-0 1 < -.0u2 Oft > 

AVERAGE PHASE SHIFT (RADIA'^S) s 9. 33SSF54E-d2 
optical transmission = 9.8a729lCE-Ql f -.C06 OB I 

FOCUSING EFFICIENCY * 6, 2 294 4 9 IE-01 < -.646 OP I 

NON-HETEROOVllf DETECTION i.-FICIENCY = 6. 2169905E-J1 C -.852 OB 

L.O. ILLUMINATION EFFICIENCY s 6. 7664677E-01 ( -.562 OR > 

maximum ANTENNA GAIK : 1. J 2 1 4 1 19E ♦ 0 9 ( 93.C92 OB ) 

RECEIVER EFFICIENCY TO I.F. = 6. 2 973106E-C1 i -2.306 OB i 
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NULTIPL" OtTECTORS 
SINP 

NOET s <», 

lOETECd.l) : 1>0*C 
IOETEC(l,2) s 2.Q«C 
IOETEC(l>3) : 
I3£TEC(1»U) s 2»0,0 
SEND 


.OtOf 
.OfOt 
tlOf j» 
1 1 Q « u t 



HULTIPL-: OiTfCTORS 


0£T£CT0O CJMPUTaTIONS, 


OdTECTOR NO. = I 

RECEIVED SlGMftL RtAi IS SHIPTEO » = 0. 

WITH RESPECT TO CENTER OF DETECTOR. 

LOCAL OSCILLATOR =iEAM IS SHIFTED X = C. 

with respect to center of DETECTOR. 

CIRCULAR JiTfCTOR DIAMETER e 1. 2962wfll 

signal POHE3 = 8.iieS9-.aZ-C7 
local oscillator POMER = I.283a6q5£-Q5 

CROSS FROruCT POWER = 2. 6^25 16 hE- Jd 

PHASE MATCH EFFICIENCY s 9. 962155 SE-0 1 ( -.002 09 > 

AVERAf.E PHASE SHIFT (RADIA95J = 9. 3388fc54E-02 

OPTICAL TRANSMISSION = 9 . 9872910E-Q 1 ( -.t06 09 I 

FOCUSING efficiency = 8. 229449 JE-Jl ( -.846 OB ) 

NON-hETEROOTNE detection efficiency = 8.2189905E-C1 < -.852 09 » 

L.O. illumination EFFICIENSY s 6.7864677E-01 ( -.562 09 ) 

maximum antenna gain = 1. ;214119E»09 < 90.C92 OB I 

RECEIVER EFFICIENCY TO I.F. = 6. 2973196E -Ci i -2.008 OB > 


AND T = C. 
ANO Y s 0. 


DETECTOR NO. * 2 

RECEI-VEO SIGNAL BEAM IS SHIFTED X s Q. ANO Y s 0. 

WITH RESPECT TO CENTER OF DETECTOR. 

LOCAL oscillator 9f AM IS 3HIFTEO X = C. AND Y s 0. 

WITH RESPECT TO CENTER OF 0£TECT0R. 

RECTANGULAR DETECTOR OIM-NSIONS ARE X * 1.2962w614E-dl AMD Y = 1 , 296 2**8 14£ -01 

SIGNAL POWER s 8.29266631-07 
LOCAL OSCILLATOR POWER * 1 . 7 8C 1969E-05 

CROSS PRODUCT POWER « 2. 894d692e-06 

PHASE MATCH EFFICIENCY a 9. 3992926E-0 1 < -.027 09 > 

AVERAGE PHASE SHIFT (RAOIAySt » 5 • 433*»2S 7E-02 

OPTICAL TRANSMISSION s 9. 987291JE-0 1 < -«C06 OB ) 

FOCUSING EFFICIENCY « 8.3J32l89E-n < -.8C8 06 ) 

NON-METEROOYNE DETECTION EFFICIENCY * 8. 2926663E-0 1 ( -.$13 09 ) 

L.O. ILLUMINATION EFFICIENCY a 7 . 90520 18E-01 < -1.021 09 ) 

MAXIMUM ANTENNA 6AU a 1. j2l4119E«>09 ( 90*092 08 ) 

RECEIVER EFFICIENCY TO I.F. a 4. 5783655E-0 1 ( -3.393 OB ) 


DETECTOR NO. * 3 

RECEIVED signal BEAM IS SHIFTED X a -6.48124L70E-02 ANO V s 0. 
WITH RESPECT TO CENTER OF DETECTOR. 

LOCAL OSCILLATOR 9EAM IS SHIFTED X a -6. 48124070E-02 ANO Y a 0. 
WITH RESPECT TO CENTER OF OETECTOR. 

CIRCULAR OETECTOR DIAMETER a 1. 29624814E-01 

SIGNAL POWER a 3.E7283U1E-07 
LOCAL OSCILLATOR POWER a 1 . 268151 8E-0 5 
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CROSS HRIJUUd HOWtK = ^ 

OH4SC match ►FFICier.'CV = i.992Jfll^£-01 < -.^^2 09 t 

AVERAGE OHAS? SHIFT (RAOIflVSI = -2.6T84753E«02 
optical TFANSKISSIOK = 9. 9672910E-01 < ••CQ6 09 I 

FOCUSING efficiency s 3. b /750 39E- J I ( -4.144 06 > 

NON-HETE?ObYM£ OcTECTtON -“FICIEHCY = 3. 6724301c-dl C -4.S5C 09 » 

L.O. ILLUMINATION EFFICIENCY a /.U73252£-0t t -1.477 06 ) 
MAXIMUM, antenna GAIN a L . 0 2 14 i 19E»0 9 ( 9i).C92 06 > 

RECEIVER EFFICIENCY TO I.F. = 6. 68087 47E-G2 ( -11.752 06 1 


DETECTOR NO. = *. 

RECEIVED SIGNAL BEAM IS SMIFTEl t * -6.48l24C7QE-tt2 AND Y a 0. 

WITH RESPECT TO CENTER OF OETECTQR. 

LOCAL oscillator REAM IS SHIFTED X = -t . 4812407 j£-02 ANO Y a Q. 

WITH RESPECT TO CENTER 0^ DETECTOR. 

rectangular OETECTCR dimensions are X a 1.29624614E-IJ1 AND V a l. 296 24814E-01 

SIGNAL POWER a 4 . 7o 1 233b E- u 7 
LOCAL oscillator POWER a l. 7595453E-03 
CROSS PPOCUCT power = 1.25jC753£-36 

phase match EFFICIENCY a b . 29l 2u3 lc-0 I i -.201 09 I 
AVERAGE OHASE SHIFT (RAQIAXSI = - 1. 68876u lt-93 
OPTICAL TRANSMISSION a 9 . 9ft 7 29ld E-O 1 < -.C06 OB • 

FOCUSING EFFICIENCY a H.7b7292J£-n < -3.217 06 I 

NON-HETEROOYNE OETECTION iFFICIENCY = 4 . 76123 35E-0 1 C -3.22! 09 J 

L.O. illumination efficiency a 6.3650 314E-01 ( -1.634 09 I 

MAXIMUM ANTENNA GAIN a 1 . d 2 1 4 1 1 9E«^ 0 9 ( 90.C92 06 ) 

RECEIVER EFFICIENCY TO I.F. a 3 . 8812059E-02 < -10.515 06 » 


177 



